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During the term of support by NASA of this research program, we have 

successfully developed a new approach to determining the radiative 

properties of complex cloud fields. This innovative research began as a 
unique, untried concept and has matured into a verified, useful 

technique. The accomplishments of this research are highlighted below. 
Following the highlights are brief explanations of each accomplishment. 

Accomplishments : 

o Design and construction of a laboratory apparatus (CFOS) to 
simulate the interaction of cloud fields with visible 

radiation. 

o Verification of CFOS by comparing experimental results from it 
with calculations performed with a Monte Carlo radiative 
transfer model. 

o Development of a software library to process, reduce and 

display CFOS data for use in research studies. 

o Utilization of CFOS to study the reflected radiance patterns 

from simulated cloud fields. 

1.0 Design and construction of a laboratory apparatus (CFOS) to simulate the 
interaction of cloud fields with visible radiation. 

The CFOS apparatus is designed to measure the simulated reflection 
of shortwave radiation by optically thick clouds. The CFOS basically 

consists of a light source, a target cloud field, an array of radiation 
detectors and a data display/collection system. Fig. 1 shows a schematic 
of the CFOS which is housed in a 10 m X 10 m dark room at the Department 
of Atmospheric Science at Colorado State University. Each of the 
components which comprise the CFOS has been selected or designed to 
simulate as closely as possible the actual interaction of visible solar 
radiation with real clouds and to perform measurements on the simulation 
analogous to the real world case. Each aspect of the CFOS is described 
below along the rationale used in selecting various options in the 
components design. 

A. The Source 


In simulating the interaction between the visible solar irradiance 
and clouds it is desirable to duplicate as closely as possible three 
characteristics of the natural source. First, the source should possess, 
to the extent possible, a visible spectrum of a 5700 K blackbody 
radiator. Second, the beam irradiance should be uniform across the area 
of the cloud field, and third, the beam should be parallel. Because the 
cost of constructing a solar simulator approaching the above description 
would have been prohibitive, efforts were focused on achieving the latter 
two characteristics for the following reasons. First, theoretically, the 
spectral distribution of the incident radiation is important only in 
relation to the microphysical properties of the cloud; i.e. scattering 
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Figure 1 . A schematic drawing of the Cloud Field Optical 
Simulator (CFOS) designed to simulate the 
interaction of visible solar radiation with optically 
thick clouds. 
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properties at a wavelength A are determined only as a function of the 
ratio A/r, and where r is the radius of the droplet. Since simulated 
clouds used in the CFOS are not composed of droplets (much less the 
specific droplet distributions measured in real clouds) and in view of 
the considerations given above, achieving the proper spectral response 
was given a lower priority and emphasis was given to designing a source 
which fulfilled the remaining criteria. Maintaining the proper spatial 
characteristics of the beam, as specified in the second and third 
criteria listed above, assures that the effects of cloud geometry are 
accurately taken into account. Doing so, even at the expense of 
achieving less than optimum spectral characteristics, is consistent with 
recent findings concerning the importance of cloud geometry. To these 
ends it was found that an array of low voltage display lamps provided 
illumination with the best spatial uniformity over the target cloud 
field. The low voltage design of the lamps allows small filament 
assemblies to be precisely placed at the focal points of accurately 
molded reflectors. The result is a smooth illumination curve which 
decreases slowly ( 7 % in the first 4°) as a function of the angle from the 
beam axis of each lamp. Use of an array of overlapping beams, each aimed 
at a different section of the cloud field, results in an irradiance 
uniform to within ± 3 % of the mean value in the plane of the cloud field. 
The source was placed behind a circular aperture 0.5 m in diameter to 
shield the cloud field from stray light and to constrain the beam to 
within 3° of its axis, thus assuring nearly parallel light. The beam is 
reflected by a large plane mirror in order to allow enough distance 
(20 m) for the nearly parallel beam to diverge to a cross section 
sufficient to illuminate the entire target area ( ~3 m 2 ). Although the 
source output is reasonably constant in “time, it is monitored by a wide 
f ield-of-view detector near the cloud field so that the effects due to 
even small variations may be accounted for in the analysis. 

3. The cloud field assembly 


The target cloud field assembly consists of the simulated clouds and 
the hardware necessary to support and orient the clouds with respect to 
the incoming radiation and the detectors. The assembly is contained in a 
separate 4mX4mX4m black box in the laboratory room. The types of 
materials used to construct the simulated clouds are discussed in 
Sections 3a-c. As shown in Fig. 1 , the target cloud field is mounted in 
a circular area -2 m in diameter. The incoming radiation enters through 

a large aperture in one face of the box. The target cloud field is 

oriented in a vertical plane and has two rotational degrees of freedom. 

The first of these allows rotation about the vertical axis of the cloud 
field and simulates changing solar zenith angle. The second-rotation is 
about a horizontal axis through the center of and perpendicular to the 
plane of the cloud field which allows simulation of a change in solar j 

azimuth angle. Each of these rotations is driven by a separate motor and 
may be controlled at the data display station. The rotations are 
monitored to the nearest degree through the use of optical encoders and 
associated logic circuitry. Any sun-cloud geometry of interest may be 

easily simulated in the CFOS as a result of this design. 



C. Sensor array 


In order to examine the relationship between the summative and 
component properties of radiation reflected from the simulated clouds it 
is necessary to sample the radiation field from a number of angular 
viewing coordinates. The placement of the CFOS sensors with respect to 
the cloud field may be seen in Fig. 1 . Fifteen sensors are mounted on a 
semi-circular ring whose ends are attached to the cloud field apparatus 
at points above and below the cloud field on its vertical diameter which 
serves as an axis of rotation for the ring. The sensors are spaced at 
increments of 10° of arc along the ring, beginning at the center of the 
ring and proceeding toward its ends. This sensor configuration permits a 
density of angular sampling of the reflected radiation which is more than 
adequate to establish the relationships between the desired quantities. 

Each sensor consists of a high-quality silicon photodiode with a 
built-in high-gain operational amplifier circuit. The diodes are mounted 
in small cylindrical enclosures each of which contains a spectral filter 
and two circular diaphragms. The diaphragms limit the field of view of 
each photodiode to 7.5° half-angle. The relative spectral responses of 
the detectors with and without the optical filters are shown in Fig. 2. 
The filtered spectral response corresponds closely to that of the human 
eye. Thus, features visually observed should be evident also in the 
measurements. The field of view of the detectors is such that the center 
detector resolves a circular area 0.5 m in diameter when viewing the 
cloud field at the nadir. By using the relationships between the optical 
and geometric properties of simulated and real clouds (see Sections 3a-c) 
it is possible to associate the CFOS sensor configuration with that of a 
real world sensor 100 km above the earth's surface with a nadir ground 
resolution of 25 km. 

D. Data display/collection station 


One of the advantages of a laboratory investigation of the radiative 
properties of clouds lies in the ability to observe an unchanging cloud 
field from a variety of angular coordinates. This feature is fully 
appreciated by observing the changes in the various detector voltages as 
the "sun-cloud-sensor" geometry changes. The output of each sensor 
including the source monitor may be read from digital voltmeters at the 
data display /collection station. An array of light-emitting diode bar 
graph displays indicates the voltage outputs of all 15 diodes 
simultaneously, giving an instantaneous sample of the radiance patterns. 
The simulated solar zenith and azimuth angles and the angle of rotation 
of the detector ring are controlled from and monitored at the display 
station. 

In order to examine the structure of the radiance patterns and to 
detect changes in the field or image as the clouds' geometric properties 
are changed it is necessary to assemble all measurements into a single 
representation of the reflected radiance pattern. This is accomplished 
in two stages by the CFOS. First, the data are recorded using an Apple 
II microcomputer as a smart data collection device. Various parameters 
such as sampling rates and the numbers of measurements included in 
average readings may be programmed into the data collection process. 
Second, the data from the floppy diskette in the Apple II are transferred 
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Figure 2. The relative spectral responses of the detectors used 
in the CFOS with and without optical filters. 



Figure 3« A comparison of the relative radiances reflected into 

the principal plane by cotton-simulated clouds (measured 
by the CFOS) to those calculated using a Monte Carlo model 
for water clouds. Comparisons are shown for solar zenith 
angles of 60° and 0 °. (Only one of the symmetric branches 
is shown in the 0° case) . Note that the uncertainty in 
observed and modeled relative radiances are a function of 
the observation zenith angle. The error brackets shown are 
maximum values typical at y = 0.55. Uncertainties near 
y — 1.0 are about half these maximum values and are indicated 
by the error bars. 



to an Apple III microcomputer where they are processed. The processing 
is relatively simple; measurements are rotated into a nadir-relative 
azimuth angle coordinate system, voltages are divided by relative 
sensitivities, effects caused by variations in, the source are removed, 
the reflected radiances are integrated over the lower hemisphere to 
obtain a flux density, and various graphical displays are generated. The 
total time of the processing of the data is -15 min. Thus, experiments 
may be repeated for slight perturbations in the cloud field geometry and 
changes in the radiance field may be appreciated within a short time. 
Also, data and displays from the radiance patterns may be saved on 
diskette files for later viewing and comparison. 

The CFOS apparatus described above should provide insight into the 
composition of wide field of view satellite measurements and images. It 
offers an interactive setting in which cloudy scenes may be created, 
observed visually, and whose relative radiance fields may be measured 
quantitatively. The next section describes the verification of the CFOS 
and in doing so furthers the appreciation of the similitude between the 
CFOS and the real world. 


2.0 Verification of CFOS by comparing experimental results from it with 
calculations performed with a Monte Carlo radiative transfer model. 

The previous sections have described a laboratory device designed to 
measure the simulated radiative properties of clouds of a realistic 
shape. In the sections which follow the properties of the simulated 
clouds and the relationship of their scattering properties to those of 
simulated water clouds are discussed. The verification of the CFOS is 
based not on analogies between the microphysical properties of real and 
simulated clouds, but rather on comparisons of observable radiative 
properties. 

A. Selection of suitable materials for cloud simulation 


The validity of measurements obtained from the CFOS depends, almost 
entirely, on the measurable radiative properties of the simulated clouds. 
Several materials were examined in a search for one which adequately 
simulates the visible reflective properties of optically thick clouds. 
The criteria which were examined in the selection process were first, the 
behavior of the radiances reflected into the principal plane by a 
horizontally semi-infinite cloud, and second, the visual appearance of 
the simulated clouds and cloud fields. Based on these criteria, two 
materials emerged as superior to all others tested; they are surgical 
quality sterilized cotton and decorative billet (DB) styrofoam which is 
marketed by Dow Chemical under product code #81568. 

It should be emphasized that the manner in which light interacts on 
the microphysical scale with the materials selected may differ 
considerably from the classical picture of independent Mie scattering in 
a spherical polydispersion. For example, the surgical cotton used in the 
experiments consists of long fibers of mainly cellulose from 10 to 25 ym 
in diameter. The mean index of refraction of cellulose is 1.55 in the 
mid-visible. The DB styrofoam is a polystyrene extrusion. It contains 
no pigments or other suspensions and has a mean cell dimension of 1 .9 mm. 
The index of refraction of polystyrene is 1.59-1.60 in the visible. 
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Figure 5. Photographs of a styrofoam-simulated cloud system. 

The bottom photograph is an enlargement of the first 
showing the detail of the area between the simulated 
cumulonimbus clouds. 


10 


Polystyrene is a transparent plastic. The white appearance of styrofoam 
results 'from multiple reflections of the incident light from the many 
facets which surround the individual cells. Thus, it would be difficult 
to draw an analogy between the local scattering process in water clouds 
and the multiple reflection process in styrofoam. Moreover, recent 
studies have shown that cloud shape may dominate microphysical scattering 
as the most important factor affecting the radiance pattern. Also, it is 
doubtful that any material which closely duplicates the scattering 
process in water clouds would possess an extinction coefficient great 
enough to permit simulation of optically thick clouds in a laboratory 
setting. The goal of the CFOS experiment is to determine if it is 
possible to simulate the observable properties of optically thick clouds 
of a realistic shape and not necessarily to duplicate the exact 
scattering process on the microphysical scale. Thus the determination of 
which materials are best qualified as simulation materials is based on 
comparison of observable optical properties of the simulation materials 
with those predicted by radiative transfer models of real water clouds. 
Such comparisons have been made for "horizontally infinite" and cubic 
finite cloud shapes. 

Figures 3 and 4 show comparisons between the radiances reflected 
into the principal plane by horizontally infinite cotton and styrofoam 
simulated clouds compared with theoretical curves derived from a Monte 
Carlo radiative transfer calculation. The model was run for a C. 1 
( Deirmendjian, 1969) droplet distribution with an absorption for a 
wavelength of 0.7 ym and for solar zenith angles with cosines of u Q = 1.0 

and u = 0.5. The radiances have been normalized to the value at the 
zenith and are displayed as a function of the cosine of the zenith angle 
(y) and the relative azimuth angle <j> into which they are reflected. 
Scattering into the forward direction (0° ^ <}> SI 5°) is depicted in the 
right halves of the figures and backscatter (165° S <p S 180°) is 
indicated in the left halves. The materials are almost equivalent in 
satisfying the first criterion. The styrofoam-simulated clouds display a 
slightly more realistic behavior in the backscatter direction. 

Both materials do well in simulating the visual appearance of real 
clouds. Fig. 5 shows photographs of styrofoam-simulated clouds of a 
realistic shape. If we add a third criterion, the ease of working with 
the materials, styrofoam is preferable for the following reasons. First, 
styrofoam is rigid and will maintain a constant ratio of geometric to 
optical depth, and the shape of simulated clouds will be retained 
regardless of orientation. Second, the simulated styrofoam clouds are 
easily placed and maneuvered in a cloud field and may be categorized and 
stored for use in other cloud scenes. Finally styrofoam clouds are 
easily made into regular shapes i.e., cubes, cylinders, spheres, etc., so 
that their radiation properties may be compared with theory. 

B . The optical depth of simulated clouds 


In order to relate the reflective properties of simulated clouds to 
real water clouds it is necessary to have an accurate estimate of the 
optical depth of the simulated clouds. One means of doing so is based on 
the bulk radiative properties of a 'horizontally-inf inite ' sheet of the 
cloud material at 0° zenith. Specifically, theory predicts that the 
ratio of the spectral reflectance to the spectral transmittance Tj of 
a semi -inf inite cloud over a non-reflecting surface is nearly a linear 
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function of optical depth (see, e.g. , Coakley and Chylek, 1975; Stephens, 
1978). Fig. 6 shows a plot of the ratio of R/T generated from the 
parameterization of Stephens (1978) for water clouds in the visible 
portion of the spectrum. The solid curve shown in the figure is the 
result of a linear regression of the form R/T = xdnx + 8.307) /1 53 • 1 49 , 

which is nearly linear for 10 < t < 200. Also shown is the ratio ( R ^/T^) 
calculated using an Eddington model for various microphysical 
distributions and a few results of the R A /T A ratio from Monte Carlo 
calculations. The nearly linear dependence between R/T and t, the 
visible optical depth, is apparent. 

An equivalent visible optical depth (t q ) may be assigned to the CFOS 
cloud material by measurement of the ratio of R/T for a sheet of the 
material with a large ratio of horizontal to vertical dimension. Once t 
is established for a large slab of styrofoam several clouds may be 
"sculptured" from the slab while maintaining in each the vertical 
dimension of the original slab. The procedure requires an additional 
stop in the case of cotton since the vertical structure of the original 
"slab" (surgical cotton is available in 12" X 60" sheets) cannot be 
maintained while forming clouds of a realistic shape. However, it was 
found that the mass of the material in a vertical column of unit cross 
section is also related to the ratio of the R/T in a nearly linear manner 
(see Fig. 7). Thus, for a cloud simulated from surgical cotton on 
approximate visible optical depth may be assigned by weighing the cloud 
and measuring its cross-sectional area to determine the area mass 
density, then using the R/T line as a transfer function to obtain t . 

C. Retrieval of basic finite cloud features 


A crucial test of the CFOS is the ability to measure the reflected 
radiances of finite clouds. Figs. 8 and 9 show comparisons between 
radiances reflected into the principal plane by finite clouds as measured 
on the CFOS and the same quantities predicted by the theoretical Monte 
Carlo model described in McKee and Cox (1976). The radiances in each 
case have been normalized by the radiance measured or calculated at the 
zenith. The CFOS profiles have been retrieved from measurements over a 
field of simulated (styrofoam) cubic clouds cut from a slab of material 
whose R/T ratio corresponded to an equivalent vertical optical depth 

T e = 76. The simulation clouds were placed at the centers of adjacent 
squares of sides 5w, where w is the dimension of the finite cubic cloud. 
This arrangement resulted in a true fractional cloud cover f =0.04. The 
cloud rows were aligned along two mutually perpendicular axes, one of 
which formed the line of intersection of the principal plane with the 
surface. Figs. 8 and 9 show measured retrievals. The first (applicable 
only for observation of the cloud field in the direction of the cloud 
rows) assumes that as the observation zenith angle (0) increases, the 
solid angle of the observed cloud field relative to the underlying 
surface increases by a factor of 1 + tan0, resulting in a modified 
fractional cloud cover f' given by 


f' 


f(l + tan0), 

f = f ( 1 + tan0 ) , 
max max 


0 < 0 
0 > 0 


max 

max 


In the above, F ^ is the fractional cloud cover derived as the tops and 
sides of the cubic clouds maximally obscure the inter-cloud spacing. 
Note, in the present case, for observations taken in the principal plane 
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Figure 6 . The ratio of reflectance to transmittance for horizontally 
infinite clouds plotted as a function of optical depth in 
the visible region of the spectrum for a solar zenith angle 
of 0°. Values of the ratio from the Eddington and Monte 
Carlo models are valid at 0.7 pm, values from Stephens (1978) 
are representative in the visible (A < 0.75 pm). 
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Figure 7 . The ratio of reflectance to transmittance for cotton 
simulated "horizontally infinite" clouds plotted as a 
function of the area mass density of the material which 
is used as a transfer function to obtain the equivalent 
optical depth of simulated clouds. 
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Figure 8- A comparison between calculated and simulated radiances reflected 
from finite cubic clouds into the principal plane for a solar 
zenith angle of 0° using styrofoam as the cloud simulation 
material. (Only one of the symmetric btanches is shown for the 
modeled results.) Note that the uncertainty in observed and 
modeled relative radiances are a function of the observation 
zenith angle. The error brackets shown are maximum values typical 
at jj = 0.55. Uncertainties near y = 1.0 are about half these 
maximum' values and are indicated by the error bars. 
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Figure 9.^ As in Fig. 8, but for a solar zenith angle of 60°. 
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(parallel to the row directions), f = 0 .20 and e max = 76°, which is 
larger than the zenith angle of any of the observations depicted below, 
so that the first line in the equation above applies. The second 
treatment utilized the exact geometry of the cloud field and the CFOS 
detectors, counted the clouds in the field of view, and calculated the 
total solid angle subtended by the clouds and their shadows. In this 
interpretaiton the radiance measurement (N) is given as 

cld cld w cir cir u sdw» 

where the subscripts cld, clr, and sdw represent cloud, clear and shadow, 
respectively, and Q is the fractional solid angle subtended by each 
target component. It was assumed that there was no radiance contribution 
from the shadow regions. The clear measurements were made of the surface 
on which the clouds were mounted in the absence of the clouds. The 
albedo of the surface was 0.02. The second treatment resulted in a 
"noisier" curve because a single cloud was either considered completely 
within or exterior to the field of view. Nevertheless, smoothed versions 
of the curves are in good agreement. The apparent disagreement between 
the curves for the 0° solar zenith case at large observation zenith 
angles is principally from two factors. First, the styrofoam material is 
least faithful in its reproduction of the reflected radiance of real 
clouds at the larger observation angles as evidenced in the comparisons 
for horizontally infinite clouds in Fig. 4. Second, in the Monte Carlo 
model, for a 0° solar zenith angle, there is no radiation incident on the 
vertical side walls of the cubic clouds. However, in the laboratory, the 
incident radiation is not parallel in the strictest sense. Thus, over 
the extended cloud field, some incident radiation falls on the vertical 
walls of some of the simulated clouds. This results in relatively 
brighter side walls in the simulated case which contribute greatly to the 
radiance field at larger observation zeniths. 

It should be noted that the retrieval of the finite cloud reflected 
radiances described above represents a rigorous test of the CFOS concept. 
The low surface albedo (0.02) combined with the extremely small 
fractional cloud cover (0.04) resulted in reflected radiances at the 
lowest bound of the domain for which the CFOS was designed. Larger 
clouds, smaller cloud spacing and more realistic simulated surfaces will 
all tend to increase the radiance signal. In addition, the above 
experiment represents the retrieval radiances reflected from only part of 
the scene. Measurements of total scene radiances are of a more stable 
nature. 


3.0 Development of a software library to process, reduce and display CFOS 
data for use in research studies. 

Introduction 


This manual was written in order to provide a transfer of established 
methodology used in the operation of the Cloud Field Optical Simulator 
(CFOS) to interested users of the device. Throughout this document it is 
assumed that the reader is familiar with the basic design and purpose of 
the CFOS. If this is not the case the reader is referred to a publication 
Davis, Cox and McKee (1983), which will provide sufficient backround 
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information and which may be found in Appendix A, This manual is divided 
into two parts. Part one is a description of the methodology of data 
collection. In part one the reader will find a discussion of alignment of 
the CFOS, calibration of the sensors, and the normal procedure which has 
been developed to collect meaningful data. Part two is a discussion of 
data processing. The main thrust of part two focuses on how to proceed 
from the raw data to an array of bi-directional reflectance data. This, 
of course, is not the only product for which the CFOS was designed but 
the basic methods are demonstrated in obtaining this type of data. Part 
two necessarily discusses some of the application software which has been 
written for the CFOS although additional software exists for particular 
tasks. All the software is listed in Appendix B along with a brief 
statement of its purpose. Finally, it must be noted that all of 
procedures associated with the CFOS have been in a constant state of 
flux. CFOS itself has always been in part the object of the ongoing 
research effort. Thus, especially in using the applications software, 
some user familiarity with the code must be established. CFOS is not now, 
nor in ail probability will ever become a black-box laboratory device. 

I. PART I 

A. Laboratory procedures 

Part I is divided into two sections. Section one presents 
suggestions, advice and procedures to ensure that the basic laboratory 
device is aligned for data collection. Certain care must be taken before 
any data is collected especially if there has been a significant down 
period or if it is discovered that the device has been moved or altered 
as may happen for instance if items not related to CFOS have been moved 
in or out of the basement annex. Section two discusses the calibration 
and method of data collection which has been found to provide the best 
results to date. It should be stressed at the outset that calibration 
should be considered part of the data collection process since it has 
been found that the calibration coefficients are not constants. In all of 
part one it is assumed that the reader has been briefed on the proper way 
to power up the laboratory and on the functions of the various controls 
accessed at the data display station. 

B. Section 1 


1 . Source alignment 

The source may need to be realigned periodically because of the 
reasons mentioned above. The source alignment consists of adjustment of 
the relative position and orientation of the light source itself and 
the mirror and establishing the zero degree position of the cloud field 
ring all of which are interdependent. An iterative procedure has been 
developed which seems to work reasonably well; however, it is somewhat 
time consuming. One begins by ensuring that the normal to the cloud field 
ring is pointing in the approximate direction of the light source. Then 
by adjusting the position of the array of lamps so that its lateral 
positioning behind the aperture appears symmetric when viewed by an 
observer standing on either side of the cloud field ring. In other words, 
if you can see all the lamps except for part of one of the lower lamps 
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when standing on one side of the cloud field ring, then if you move to 
the other side of the ring you should be able to see all of the lamps 
except for a similar portion of the other lower lamp. The position and 
orientation of the mirror may also need adjustment in order to accomplish 
this. It may take a few to several adjustments obtain the desired 
symmetric alignment. 

Next, the true zero zenith angle must be set. To do so it is first 
necessary to establish the center of the cloud field with respect to the 
direction of the incoming beam. An aluminium bar has been constructed for 
this purpose. It is fitted with cylindrical feet on each end which may be 
inserted into the holes in the cloud field ring. A small hole has been 
drilled in the middle of the bar through which the source light may pass 
to form a small image of the lamp source on the cloud field backround. 
The image marks the center of the cloud field with respect to the 
incoming beam. Next, an alignment photodiode which is mounted in a 
circular magnetic support is placed on the cloudfield backround at the 
center point. The photodiode is connected to an operational amplifier 
circuit and thus requires a supply votage of +/- 15 volts. The output of 
the device may be read using a digital voltmeter. By examining the output 
of the device for changes in the simulated solar zenith angle about the 
initial orientation, one may establish the center of the response 
pattern. This point marks the true zero of the solar zenith angle 
regardless of the reading on the data display panel meter. The zero 
readout may be reset once the optical zero is established. It should be 
pointed out that near the zero of the solar zenith angle there will be 
little change in the output of the device due to the nature of the cosine 
response. Thus one would be well advised to look for symmetric points at 
larger angles about the zero in solar zenith. 

After the optical zero zenith angle is established the next step in the 
general alignment procedure is to check the spatial homogeneity of the 
incident irradiance field. The magnetically held diode may be used to 
check the homogeneity of the incident radiance field. By moving the 
device about to different positions on the cloud field and comparing the 
output voltage at the various points homogeneity may be established. If 
the incident pattern is not homogeneous then the orientations of the 
individual lamps must be must be adjusted relative to each other. This 
may be a time consuming process depending on how non-homogeneous the 
field is and on how large the region is where homogeneity is desired. It 
is convenient to employ a second digital voltmeter to monitor the output 
of the test diode. The second dvm should be placed near the lamp source 
so that it can be read at the source. Then while adjusting the individual 
lamps with the test diode placed in a region of particularly high or low 
irradiance the effect of the re-aiming of the individual lamp may be 
monitored. The dvm located near the interior of the CFOS is convenient to 
monitor changes in the test diode output as the diode is moved about on 
the cloud field. Two items should be noted while checking the homogeneity 
of the incident beam. First, the test diode circuit outputs a dark 
current voltage which must be subtracted from the readings before 
establishing the relative magnitude of the incident irradiance. Second, 
the reorientation of the individual lamps changes to some extent the 
optical zero zenith angle. If in the process of smoothing the incident 
light one upsets the optical zero the preceding step should be repeated. 

The next item in this section concerns the cleanliness of the glass 
surfaces in the CFOS. The large mirror will need cleaning from time to 
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time and the frequency of the cleaning will be obvious. Also, the 
photodiode glass covers and the glass optical filters housed within each 
collimating can should be cleaned periodically, especially after 
significant down time. When, this is done the individual diodes invariably 
become misaligned and uncalibrated. Recalibration will be covered in 
section 2 of part 1. To re-align each diode a simple flashlight device 
has been construced which may be mounted over the exterior of each diode 
can and which projects a small spot of light toward the center of the 
cloud field backround. The device is powered by any 4.5 volt d.c. source. 
It is helpful to mount a white piece of paper or cardboard at the center 
of the cloud backround and to turn off the main lights since the 
projected spot of light is rather dim. The center of the cloud field may 
be established on the white reflector as explained above using the large 
aluminium bar. In this manner each diode's field of view may be centered 
at the center of the cloud field after the glass surfaces have been 
cleaned. 

The final item to be discussed in this section concerns the alignment of 
the detector arc. This is one of the simplest of the procedures. The 
optical zero of the detector ring is established at the position where 
the shadow of the detector ring falls on the center of the cloud field. 
Actually, there are five shadows since there are five lamps, but it is 
fairly oovious where one should place the detector arc zero point. The 
readout may be reset to zero at this point if desired but this is not 
necessary if the offset is programmed into the software. 

C. Section 2 


This section presents a discussion of procedures used more on a 
daily basis or even more often to ensure quality data. These methods 
pertain to the actual process of collecting the data and maintaining 
calibration. It is easiest to begin with the data collection device--the 
Apple II plus computer and the data logging program. After turning the 
computer on, with the disk labeled "cfos prodos" inserted, the logging 
program is initiated by typing "run calib." This program will prompt the 
user to insert a properly formatted data disk and to select the number of 
data scans to be included in each measurement. In an actual data 
collection session 30 data scans have usually been used. However, for 
merely observing the diode outputs a smaller number is more convenient. 
The program then displays the output of each of the diodes and the 
relevant angle settings. Nothing else will happen until the user presses 
a "control p" which initiates a write to the disk at the end of the 
current scanning sequence. Data may be taken in thi3 manner for as many 
different geometries as desired. When it is desired to end the recording 
session the user must press the "escape" key which will terminate the 
file. This is the basic way used to collect the data. Small adjustments 
to the basic process and specific cautions will be found below as they 
pertain to collection of bi-directional reflectance data, calibration 
etc. 


1. Calibration 


In all the applications for which CFOS was constructed only the relative 
sensitivities need be determined. Even so, it has been found that 
constant checks must be made on the diode relative sensitivities since - 
they may change significantly even with the temperature excursions 
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induced by the heating system in the basement annex. Thus what follows 
may seem like overkill, but it has been found that only by continually 
updating the sensitivities can one expect to obtain results valid to 
within a few percent. 

Each diode circuit has a very high amplification factor which causes 
a significant dark current signal. It is mainly this compontent of the 
signal which is temperature sensitive. In order to "predict" the 
magnitude of the dark current signal, (which is impossible to monitor in 
an experiment with light incident on the diodes), a separate diode 
circuit identical to the others has been masked off so that no light may 
enter the diode can. Note that it is important to ensure that the back of 
the collimator can is also masked off since it has been observed that 
sufficient light may enter the can from the backside and cause erroneous 
results. The purpose of this diode is to put out only a dark current 
signal which will respond to temperature changes and which may be used to 
infer the behavior of the dark current signals in the remaining diodes. 
Simple linear regression has been found to provide predictions of the 
required dark current signals which are highly correlated to the actual 
measured valued. (Correlation coefficients greater than 0.99 are not 
uncommon.) The regression relations should be re-calculated periodically 
using the following procedure. 

It is desirable to collect the regression data over a period of time 
when the temperature extremes in the CFOS laboratory are similar to those 
encountered during the data collection session. Thus it may be desirable 
to collect the data over the weekend before a series of measurements will 
be made or if the outside ambient temperature has changed drastically to 
collect a new set overnight. By modifying the data logging program 
slightly the system may be placed in a continuous logging mode and 
serveral hundred scans may be recorded from which the regression data may 
be generated. 

The modification consists of inserting the line, 

"1516 GOSUB 600" 

Depending on the amount of time the program is left to record data, 
a very large may be generated. The software which processes this data is 
described in Part II. During the time these data are being collected it 
is necessary that the curtains on CFOS remain closed and is would be 
preferable if the main overhead lighting remain off. 

Once the regression data have been collected the relative sensitivities 
may be determined. These numbers also change and must be continually 
updated. In fact, below in the description of a typical data collection 
session note that a calibration update is performed before each scene 
measurement. To obtain the relative calibration coefficients an isotropic 
reflector is used. It has been found that ordinary polyurethane foam 
provides a surface which is nearly istropic out to observation angles of 
nearly 60 degrees when illuminated at normal incidence. It is convenient 
to use magnets to mount the foam on the cloud field backround since it 
must be mounted and removed frequently. With the foam in place the 
reflected signal from all the lamps is too large for the highly sensitive 
detectors. In order to circumvent this problem a separate power supply 
has been connected only to the center lamp of the light source. A circuit 
breaker is provided to ensure that the output of this power supply may be 
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isolated at all times from the AC supply which normally operates the 
source. 

Thus using only the center lamp one may collect data and record it 
on a disk file which may be analyzed to obtain the relative 
sensitivities. It is advisable to set the detector arc angle at about 10 
to 15 degrees in order to displace the shadow it casts as far as possible 
from the center of the foam. The software to produce the diode 
sensitivities is descibed in Part II. 

2. A typical data collection session 


To begin a typical data collection session it is assumed that the 
regression data have been collected and that the relative sensitivities 
are known fairly well. What is described here is a procedure which when 
followed will produce data with standard deviations of a few percent of 
their respective mean values. 

Although dark current data have already been collected best results 
will be obtained if this information is updated prior to each experiment. 
To do so the CFOS curtains are closed and a data file is written with all 
lights off. This file need consist of only 1 0 to 20 records each 
containing 30 or so scans. While this data is being collected the lights 
should remain off and none of the motor controls should be activated 
since the motors and switches may introduce noise into the diode signals. 
End the file with the usual escape sequence. 

Next leaving the main AC power off the center lamp of the source array is 
activated by powering up the separate DC supply. The circuit breaker 
switch at the lamp source must also be turned on. Open the curtains and 
mount the calibration foam so that it is as much in the center of the 
cloud field support as possible. Offset the detector arc angle to about 
15 degrees so that the shadow will not be in the field of the detectors. 
Write a file of about 10 scans to establish current calibration 
coefficients. End the file with the escape sequence. 

Now the system is ready for data collection. Turn off the center 
lamp supply and open the circuit breaker. Turn on the main lamp switch. 
It is best to let the source lamps heat up as their output will fall 
considerably in the first 5 minutes or so. Variation in source strength 
is accounted for in the data reduction software but a constant source is 
best. Remove the foam calibration reflector and mount the cloud field 
target. It is easiest if the cloud field is premounted on a separate 
backround so that it may be rapidly mounted and removed. Normally 
bi-directional reflectance data is collected by scanning the detector arc 
through plus and minus 60 degrees with the cloud field set at the desired 
zenith angle. It is best to confine the total target to within about 9 
inches of the center in order to avoid errors intoduced by changes in the 
solid angle subtended by cloud elements as the detectors are rotated. So 
beginning with the detector arc offset at 60 degrees from the center, on 
a new data file, record a data scan. Then proceed every 10 degrees in 
detector arc across the reflectance field. The data record which is being 
averaged while the arc is moved should be not be recorded since moving 
the arc will introduce noise into the data. Only the next record during 
which no arc movement has occurred should be recorded. After some 
practice it is possible to move the arc while the computer is writing to 
di 3 k in which case the next record may be used. If, after the "control 
p" sequence has been entered, it is desired to inhibit the write to disk, 
press any key on the Apple II keyboard and the current record will not be 
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written to disk. This happens occasionally if doors are opened or in some 
cases people may need to walk across the path of the incident beam. These 
records should be discarded. It often happens that moving the arc will 
cause some change in the zenith angle setting in which case the zenith 
angle should be adjusted to the desired value. When all the data have 
been collected for a particular scene the file should be ended with the 
escape sequence. For best results the same scene should be repeated as 
many as four times to assure the quality of the data. (This should yield 
results with standard deviations of less than 2% of the mean values.) 
Each time the experiment is repeated the dark current and sensitivities 
should be redone as explained above. 

At the end of the session the curtains should be closed, the 
computer turned off and the lamps should be turned off. However, the 
power to the main control panel should be left on. It has been found that 
the stability of the entire system is highest if the amplifier circuits 
are constantly powered up. The data disk should be removed and saved for 
processing, which is the subject of Part II below. 

II. PART II 


In part two the software which has been developed for the processing 
of bi-direcional reflectance data is described. Additional software 
exists and is listed in Appendix B. Only a mention of the particular use 
of this latter software will be given since it is unlikely that a new 
user will have need for these routines without some modification. As 
mentioned in the introduction, it is assumed that the user will become 
familiar with the software. Many will require some user interaction. All 
of the programs which are discussed below are found on the Davong system 
hard disk under volume prefix /areal or .HI or on the CFOS software 
backup diskette included with this manual. 

We begin by assuming the user has completed a session in the CFOS 
laboratory and has written a data disk containing a file for establishing 
the dark current regression coefficients, one for updating the 
coefficients, a file containing relative calibration data and one or more 
files containing data for analysis of bi-directional reflectance data. 
All of the files written by the Apple II in the CFOS laboratory have a 
file name of the type D.mmddyy. hhmm, where mm is the two digit designator 
of the month, dd designates the day of the month yy the year, hh the hour 
on a 2k hour clock and mm the minute when the file was written. For 
example, the last dark current regression data recorded is listed on the 
first page of the Appendix B as D. 082484.1 348. 

The first task is to run the regression program which will establish 
the initial regression coefficients between the masked diode and the 
remaining detectors. Run the program called "CFOSREG". It will prompt you 
for the name of the file containing the regression data. It expects to 
find the file on device .d2. After the file name is input the reqression 
coefficients will be calculated and written to a file on the Davong 
system hard disk called "/areal /cfos. cal/coeff" . As the data is processed 
the values of the slope, intercept and regression coefficients will be 
written to the screen. If the file contains a large amount of data the 
processing may take several minutes to a half hour. Once the regression 
coefficients have been calculated it is a good idea to transfer the raw 
data file to the hard disk for temporary storage in case the file 
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containing the coefficients is lost as has been done with the file 
d. 082484. 1348. 

The program above is utilized only occasionally after a new 
regression is needed for whatever reason. More often the file 
"/areal /cfos. cal/coeff" will already exist on the hard disk and these 
coefficients will need to be updated. This is the task of the program 
called "/areal /cfos. cal/cfosdark" using the data collected just prior to 
the relative sensitivity measurement described in Part I. This program 
also prompts for the name of the data input file which it expects to find 
on device .d2. The program uses the regression coefficients in 
"/areal /cfos. cal/coeff " and simply updates the offset coefficients which 
it then writes to a file called "/areal /cfos. cal/adjust" . Tnis procedure 
is normally repeated for each new bi-directional reflectance data set 
prior to running the relative calibration program described next. 

Now the program "/areal /cfos. cal/cfossens" is run to process the 
relative sensitivity measurements. The user is asked for the name of the 
data file which it expects to find on device . d2 after which the contents 
of "/areal /cfos. cal/cfoscoeff" and "/areal /cfos. cal/adjust" are used 
along with the current data file to establish sensitivities relative to 
diode #8. These values are written to a file "/areal /cfos. cal/sensiv". 

The preliminaries out of the way it is now possible to process 
meaningful bi-directional reflectance data. This is the task of the 
program "/areal /cfos. bdr/cfosbdr" which expects to find the raw data in 
the subdirectory "/areal /cfos. back/" . If several bi-directional 
reflectance files have been collected it is convenient to use the filer 
to transfer these files to the hard disk prior to the data processing 
session. The present routine uses the data in the files 
"/areal /cfos. cal/coeff ", "/areal /cfos. cal/adjust", and 
"/areal /cfos. cal/sensiv" along with the present raw data to establish the 
normalized bi-directional data set. Normally the regime of detector arc 
angles is limited to +/- 60 degrees in these data. The program terminates 
by writing to a file in the subdirectory "/areal /cfos. bdr/" concatenated 
with the current raw data file name. The bi-directional field, the number 
of detector arc angle scans, the normalizing flux density and the minimum 
and maximum of the field are also written to the file for later use. 

This completes the main processing for bi-directional reflectance 
data. However, now that the field has been established it is fairly 
common to continue the process further to obtain a denser interpolated 
field which may be plotted or used for error analysis, for example. So 
continuing along the data processing chain of events one may run the 
program "/areal /cfos. interp/cfosinterp" which performs a weighted 
interpolation on the normalized bi-directional reflectance data. The 
weighting is based on the inverse of the angular distance between the 
point at which the interpolation is desired and a neighborhood around 
this point which may be changed to include all of the normalized data 
points or just a surrounding few. The running time is highly dependent on 
the number of points used in the interpolation. The program prompts for 
the number of nadir and azimuthal points at which interpolated results 
are desired as well as the input data file name. This latter file name 
should be the concatenated portion of the normalized data which was 
written by "/areal /cfos. bdr/cfosbdr" to the subdirectory 
"/areal /cfos. bdr/" . This is of course just the "d. mmddyy. hhmm" designator 
of the original raw data. The product of this program is an interpolated 
field which is written to a file under the subdirectory 
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"/areal /cfos. plot/" concatenated with the "d. mmddyy. hhmm" designator of 
the original data. 

One final product of this line of processing is a polar plot of the 
bi-directional reflectance data. Such a plot is the product of a program 
"/areal /cfos. plot/cfosplot". This program utilizes the HP 7 470A plotter 
to produce the drawing. The input data set is the product of the previous 
interpolation program and is requested by the program. Another use of the 
interpolated field is variance estimation. It is desirable from time to 
time to analyze the variance in the bi-directional reflectance data for 
the same scene measured several times or the variance among sets of data 
for slightly different scenes. This type of analysis is the object of the 


program "/areal /cfosste" 

which 

prompts 

for 

first 

the number 

of 

interpolated 

data 

fields 

and then 

the name 

of 

each one 

. Again, it 

is 

assumed that 

the 

interpolated fields are 

to 

be found 

in under 

the 


subdirectory "/areal /cfos. plot/" with concatenations of the raw data file 
names to be input by the user. The program presently assumes that 144 
interpolated data points are available for each scene in a 12 by 12 array 
over the azimuth and nadir regime. The poduct is the standard error at 
each grid point and an average standard error for all grid points. This 
program was used to establish the basic precision of the bi-directional 
reflectance analysis which has been outlined in this manual. 

This completes the description of the software needed to obtain 
bi-directional reflectance information. There are several other programs 
listed in Appendix B. What follows is a very brief description of the 
remaining software. 

"/areal /cfos. cal/cfoscal" was the predecessor to the current 
sensitivity procedure. It assumed the dark current coefficients were 
constants and were entered in DATA statements within the program. 

"/areal /cfos/cldareadist" is a simple algorithm which implements a 
cloud field statistical model given by Planck. 

"/areal /cfos/clddistl " is a simple algorithm which implements the 
cloud field distribution program by Hozumi. 

"/areal /cfos/cfosmovie" was used to provide a "movie" type graphical 
display of bi-directional reflectance fields. The input files to this 
program were previously generated fotofiles and are assumed to exist on 
device . hi . 

"/areal /cfos/cfosf in" was used to produce graphical displays on the 
monitor and a pricipal plane graph on the silentype of bi-directional 
reflectance patterns for finite cloud fields. 

"/areal /cfos/cfosback" was originally used experimentally to adjust 
for the backround reflectance pattern but was not found to be useful. 

"/areal /cfos/cfosoff" was used to analyze the effect of displacing a 
single cloud or a small field clouds from the center of the CFOS cloud 
field support. 
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"/ areal /plot. pack/myplot", "/areal /plot. pack/myplotl " , and 
"/areal /plot. pack/hp. text " were prototype plotting routines from which 
the cfos plotting routines were patterned. 

"/areal /plot . pack/rectterp" and "/areal /plot . pack/recplot " are 
routines analagous to their cfosinterp and cfos. plot 
counterparts( described above) and were used originally to interpolate 
and plot the radiance patterns exiting modeled and measured finite cubic 
clouds. 


References 

Coakley, J. A., and P. Chy’lek, 1975: The two stream approximation in 

radiative transfer: Including the angle of the incident radiation. J. 

Atmos. Sci. , 32, 409“4l8. 

Davis, J. M. , S. K. Cox, and T. B. McKee, 1983: Design and verification of a 

cloud field optical simulator. J. Appl. Meteor., 22, 947-958. 

Deirmendj ian, D. , 1969: Electromagnetic Scattering on Spherical 

Poly-dispersions. Elsevier, 290 pp. 

Stephens, G. L., 1978: Radiation profiles in extended water clouds. II: 

Parameterization schemes. J. Atmos. Sci., 35, 2123-2132. 



34 


List of Figure Legends 

Figure I. A ohotograph of the radiance pattern from an optically thick 
cloud under normal irradiation as simulated m the CFOS. 
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6.0. Statement on Inventions 

There were no patentable inventions made during the performance of this work. 
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APPENDIX 


Data Reduction and Analysis 


Programs for the Cloud Field Optical Simulator 



A- 1 


10 REM PROGRAM CFOS 

20 PRINT" CFOS CALIBRATION PROGRAM" 

30 DIM AT (15) 

40 DIM VALUE (500,20) 

57 DIM DC (19) 

70 DIM SLOP ( 15) , OFFS ( 15) , SENS ( 15) 

75 DATA -633.57,66,195,-357,-240,-52,-59,21,-1,-84, -221,-192,-182,307 
80 DATA 1 . 3758,0. 811,0. 839,0.694,1 . 486, 1 . 049, 1 . 1 49 , 0. 78 , . 85 , 0 . 89 , 1 . 05 , 1 . 1 

.215, 1.041 ,0.597 

100 FOR K=1 TO 15 
110 READ OFFS ( I- ) 

120 NEXT h 

130 FOR H = 1 TO 15 

140 READ SLOP (k) 

155 NEXT h 

190 DTRD=3. 14159/130. 

250 OPEN# 1 AS INPUT, ".D2/D. 091784. 1247" 

251 T $= "TEST . DARK " 

260 BT= " 11 

280 INPUT#1,1;B$ 

370 NS=VAL(BT) 

380 PRINT NS 
400 FOR K=1 TO NS 
410 BT=" " 

420 I NF'UT # 1 , K+ 1 ; BT 

430 F’OS 1 = 1 

435 FOR J=1 TO 20 

440 SPACE=F'0S1 

450 F'0S2= I NSTR ( BT , " ", SPACE) 

460 N0C=F’0S2-F‘0S 1 

470 VALUE ( k , J ) =VAL (MI DT ( B$ , F'OSl , NOC) ) 

480 F'OSl =F'0S2+1 

490 PRINT VALUE (K,J) 

500 NEXT J 

510 NEXT h 

520 FOR k= 1 TO NS 
530 FOR J=1 TO 15 

540 SENS ( J ) =SENS ( J ) + ( VAL ( k , J ) *SLOP ( J ) +OFFS ( J ) ) / NS 

550 NEXT J 

560 NEXT K 

570 FOR K=1 TO 15 

580 SENS ( J ) =SENS < J ) /SENS (8) 

590 NEXT K 

600 0F'EN#4 AS OUTPUT ,". SILENTYF'E" 

610 OUTPUT #4 

620 FOR L=1 TO 15 

630 PRINT SENS (K) 

640 NEXT F. 

650 END 



10 REM PROGRAM CFOS 

20 PRINT" CFOS DARK CURRENT UPDATE PROGRAM" 

30 DIM A* ( 15 ), VALUE < 100,20) 

40 DIM SLOP ( 15) , OFFS < 15) 

50 PRINT" INPUT FILE NAME OF DARK CURRENT UPDATE RUN" 

60 INPUT D4> 

250 OPEN# 1 AS INPUT , " / AREA1 /CFOS . CAL/COEFF " 

260 FOR J=1 TO 15 

270 READ# 1 , J ; SLOP ( J ) , OFFS ( J > 

280 NEXT J 

290 0PEN#2 AS INPUT , " . D2/"+D$ 

300 8*="" 

310 I NF'LJT #2 , 1 ; NS 

320 MS=NS-2 
380 PRINT NS 

400 FOR f- = 1 TO NS 

410 B*=" " 

420 INF‘UT#2,K+1 ; B* 

430 F'OS 1 = 1 

435 FOR J = 1 TO 20 

440 SPACE=F’0S1 

450 F’0S2= I NSTR ( B-£ , " ", SPACE) 

460 N0C=F'0S2-P0S 1 

470 VALUE (K,J)=VAL<M ID* <B* , F'OSl , NOC) ) 

480 P0Sl=F’0S2+i 

490 F'F: I NT VALUE < K , 3 ) 

500 NEXT J 

510 NEXT K 

520 FOR J = 1 TO 15 

530 DIFF=0 

540 FOR K= 1 TO NS 

550 D I FF=D I FF+ ( VALUE ( K , 1 6 ) *SLOF' ( J ) +OFFS < J ) -VALUE < K , J ) ) /NS 

560 NEXT K 

570 OFFS ( J ) =OFFS < J ) -DI FF 

575 PRINT DIFF 

580 NEXT J 

600 0PEN#5 AS OUTPUT , " /AREA1 /CFOS. CAL/ ADJUST " 

610 FOR J=1 TO 15 
620 WR I TE#5 , J ; OFFS ( J ) 

630 NEXT J 

640 CLOSE 
650 END 



10 REM PROGRAM CFOS 

20 PRINT" CFOS DARK CURRECT OFFSET PROGRAM" 

30 DIM AT (15) 

40 DIM VALUE (650,20) 

50 DIM SLOP ( 15) , OFFS ( 15) 

100 PRINT" INPUT THE NAME OF THE DATA INPUT FILE" 

110 INPUT DT 

120 DTRD=3. 14159/180. 

150 OPEN# 1 AS INPUT, "/. D2/"+DT 
160 B T= " " 

180 INF'UT#1 , 1 ; BT 
190 DTRD=3. 14159/180. 

270 NS=VAL(BT) 

280 PRINT NS 
300 FOR K= 1 TO NS 
310 BT= “ " 

320 I NF'UT # 1 , H'+ 1 ; BT 

330 P0S1=1 

335 FOR J=1 TO 20 

340 SPACE=F'0S1 

350 F'0S2= INSTR ( B $ , " ", SPACE) 

360 N0C=F’QS2-F'0S 1 

370 VALUE (F, J) =VAL (MID* (BT,P0S1 ,NOC) ) 

380 P0S1 =F‘0S2+ 1 

390 PRINT VALUE (K,J) 

400 NEXT J 

410 NEXT K 

415 FOR J=1 TO 15 

420 N=0 

430 SX=0 

440 SY=0 

450 SXX=0 

460 SYY=0 

470 SXY=0 

480 FOR K= 1 TO NS 

490 N=N+ 1 

500 X = VALUE <K, 16) 

510 Y= VALUE <K , J ) 

520 SXY=SXY+X*Y 

530 SX=SX+X 

540 SY=SY+Y 

550 SXX=SXX+X*X 

560 SYY=SYY+Y*Y 

570 NEXT K 

580 BETA1= (N*SXY-SX*SY) / (N*SXX-SX*SX ) 

610 BETA0=SY/N-BETA1 *SX/N 

620 PRINT BETA1 , BET AO 

630 R= (N*SXY-SX*SY) /SDR( (N*SXX-SX*SX ) * (N*SYY-SY*SY) ) 

640 PRINT R 

650 SLOP ( J ) =BET A 1 

660 OFFS ( J ) =BETA0 

670 NEXT J 

680 0F’EN#3 AS OUTPUT ," /AREA 1 /CFOS. CAL/COEFF" 

690 FOR J=1 TO 15 

700 WR I TE#3 , J ; SLOP ( J ) ,OFFS(J) 

710 NEXT J 

720 CLOSE 
730 END 



10 REM PROGRAM CFOS 

20 PRINT" CFOS CAL I BRATON PROGRAM" 

30 DIM A-t (15), VALUE ( 1 00 , 20 ) 

40 DIM SLOP ( 1 5 ) , OFFS (15), COFFS ( 15 ) , SENS ( 15 ) 

50 PRINT "INPUT FILE NAME OF DARK CURRENT CAL I BRATON RUN" 

60 INPUT DT 

150 QF'EN#1 AS INPUT , ” /AREA 1 /CFOS. CAL/ COEFF" 

160 FOR J=1 TO 15 

170 READ# 1 , J ; SLOP ( J ) , OFFS ( J ) 

180 NEXT J 

200 0PEN#2 AS I NPUT ," /AREA 1 /CFOS . CAL/ ADJUST" 

210 FOR J=1 TO 15 

220 READ#2 , J ; COFFS ( J ) 

230 NEXT J 

240 PRINT" CURRENT DARK CURRENT OFFSET UPDATE VALUES" 

250 FOR J=1 TO 15 
260 PRINT COFFS ( J ) 

270 NEXT J 

280 CL0SE#2 

290 0PEN#2 AS INPUT , " . D2/ "+DT 
300 

310 INF'UT#2 , 1 ; NS 
320 NS=NS-2 
380 PRINT NS 
400 FOR h=l TO NS 
4 1 0 B*= " " 

420 I NPUT #2 , K+ 1 ; BT 

430 F‘05 1 = 1 

435 FOR J=1 TO 20 

440 SPACE=P0S1 

450 P0S2=INSTR <B$, " ", SPACE) 

460 N0C=F'0S2— F'OS 1 

470 VALUE <k , J ) = VAL (M I D-f (BS , F'OSl , NOC) ) 

480 F'QS1=PQS2+1 

490 PRINT VALUE <K , J ) 

500 NEXT J 

510 NEXT K 

520 FOR J=1 TO 15 
530 SENS ( J ) =0 . ' 

540 FOR K= 1 TO NS 

550 SENS (J) =SENS (J)+ (VALUE ( K, J ) -VALUE <K, 16) *SLOP<J> -COFFS (J) ) /NS 

560 NEXT K 

570 PRINT SENS ( J ) 

575 NEXT J 

580 FOR J=1 TO 7 

590 SENS < J ) =SENS ( J ) /SENS (8) 

600 SENS (J+8) =SENS ( J+8) /SENS (8) 

610 NEXT J 

620 SENS (8) =1 

625 0F'EN#5 AS OUTPUT ," /AREA1 /CFOS. CAL/SENSIV" 

630 FOR J=1 TO 15 
640 WR I TE#5 , J ; SENS ( J ) 

650 NEXT J 

660 CLOSE 

670 OPEN# 1 AS OUTPUT, ".SI LENTYPE" 

680 FOR J=1 TO 15 
690 F‘R I NT # 1 SENS ( J ) 

700 NEXT J 

710 CLOSE 
720 END 



A- 5 - 

10 PRINT" ENTER THE TOTAL CLOUD COVER" 

20 INPUT ST 

30 PR I NT "ENTER THE CHI PARAMETER" 

40 INPUT CHI 

50 PR I NT "ENTER THE DIAMETER LIMITS D1 AND D2" 

60 INPUT D 1 , D2 

70 PR I NT "ENTER THE EXPONENT ALPHA" 

80 INPUT ALPHA 
100 D 1 A=D 1 *ALPHA 

110 D2A=D2* ALPHA 

120 A=- (ST / < 2*-CH I ) * (EXP (-D2A) * < D2A 2+2*D2A+2 > -EXP (-D1A) * <D1 A'"2+2*D1 A+2) ) ) 

130 PRINT D1,D2,A 
140 END 



10 REM CLOUD DISTRIBUTION PROGRAM A LA HOZUMI 
20 DIM D ( 50) ,N<50> ,A(50> ,V(50) 

25 OPEN# 1 AS OUTPUT, ".PRINTER" 

110 FOR M=5 TO 95 STEP 10 
120 A=124*EXP <— 4. 7 *M*. 01 ) 

130 B=4. 5*EXP <-3. 5*M*. 01 ) 

135 ASUM=0 

140 FOR 1=1 TO 20 

150 D 1 = I *0 . 50—0 . 50 

160 D2= I *0 . 50 

170 N ( I ) =A/B* (EXP (-B*D1 ) —EXP (— B*D2) ) *40 

175 DB 1=D1*B 

176 DB2=D2*B 

180 A ( I ) =3. 14 159*A/ (4*B-'3> * (EXP (-DB1 ) * (DBl*DBl+2*DBl+2> -EXP (-DB2) * (DB2 

B2+2*DB2+2 > ) / 100 
190 ASUM=ASUM+A ( I ) 

280 PRINT# 1M , D 1 , D2 , N ( I ) , A ( I ) 

290 NEXT I 

295 PRINT#1ASUM: F'RINT#1 : PR I NT# 1 : 

300 NEXT M 

310 END 


1 500 

0PEN#2, ' 

. GRAF IX" 

1510 

I NVO) E" 

Dl/BGRAF. INV" 

1520 

PERFORM 

INITGRAFIX 

1530 

PERFORM 

GRAFIXM0DE<7.1 ,7.1) 

1540 

PERFORM 

F ILLF'ORT 

1550 

PERFORM 

GLOAD .".HI /FOAM . 60DEGGRAF " 

1 560 

PERFORM 

GRAF I XON 

1570 

PERFORM 

GRAF IX MODE (7.1 ,7.2) 

1580 

PERFORM 

GLOAD .".HI /FOAM . 30DEGGRAF " 

1590 

PERFORM 

GRAF I XON 

1 600 

PERFORM 

GRAF IX MODE (7. 1,7.1) 

1605 

PERFORM 

INITGRAFIX 

1610 

PERFORM 

GLOAD. " . H 1 /FOAM . ODEGGRAF " 

1620 

PERFORM 

GRAF I XON 

1625 

FOR k=l 

TO 1000: PERFORM GRAF I XON: NEXT K 

1 630 

PERFORM 

GRAF IX MODE (7.1 ,7.2) 

1640 

PERFORM 

GLOAD. " . Hl/FOAM. 60DEGPL0T" 

1650 

F'ERFOF;M 

GRAFIXON 

1660 

PERFORM 

GRAF IX MODE (7.1 ,7.2) 

1670 

PERFORM 

GLOAD .".HI /FOAM . 30DEGPL0T " 

1680 

PERFORM 

GRAFIXON 

1690 

PERFORM 

GRAF I X MODE (7.1 ,7.1 ) 

1700 

PERFORM 

GLOAD. ". Hl/FOAM. ODEGPLOT" 

1710 

PERFORM 

GRAFIXON 

1720 

FOR k= 1 

TO 1000: PERFORM GRAF I XON: NEXT K 

1 730 

GOTO 15' 

40 

1740 

END 




10 REM PROGRAM CFOS 

20 PRINT" CFOS DATA REDUCTION PROGRAM" 

30 PR I NT: PRINT: PR I NT" THE DIODE SENSITIVITIES USED IN THIS PROGRAM ARE" 

40 DIM DIANG(15) , SENS (15) ,DC(19) , NADANG ( 500 ) ,RAZANG(500) , VOLT (500) ' 

50 DIM D I O ( 1 9 ) , ANG ( 3 ) , NR ( 36 ) , NRN ( 36 ) , ME AS ( 500 , 2 ) 

55 DIM AT (15) 

60 DATA 20 , 30 , 40 , 50 , 60 , 70 , 80 , 90 , 1 00 ,110,1 20 , 1 30 ,140,150,1 60 

70 DATA 1 . , 1 . 06 , . SS , . 94 , . 84 , . 95 , 1 . 1 8 , . 79 , . 94 , 1 . 03 , . 86 , . 98 , . 98 , 1 . 06 , 1 . 04 

75 DATA 1098, 125,477,852,427,573,314,338, 193,652,185,503,73,220,287,25,3 


80 
90 
1 00 
1 10 
120 
130 
140 
143 
150 
160 
170 
ISO 
190 
200 
210 

220 
•3. 'ncr 

250 

251 


260 

270 

280 

290 

300 

310 

320 

340 

350 

360 

370 

380 

400 

410 

420 

430 

440 

450 

460 

470 

480 

490 

500 

510 

520 

530 

535 

540 

545 

550 

560 


FOR 1=1 TO 15 
READ DIANG(I) 

NEXT I 

FOR 1=1 TO 15 
READ SENS ( I ) 

NEXT I 

PRINT" DIODE NUMBER DIODE ANGLE 

NF-8 

FOR 1=1 TO 15 

PRINT USING 170; I , D I ANG < I ) , SENS ( I ) 

IMAGE 10X , 2# , 14X , 3# , 16X , 1# . 2# 

NEXT I 

DTRD=3. 14159/180. 

REM DEFINE ARCOS FUNCTION 

DEF FN ACOS ( X ) =ATN (SOR < i-X*X ) /X ) 

L=0 

AVZEN=0 

OPEN# 1 AS INPUT, D2/BACKZ60’' 

T T="FOAM" 

FOR LL=1 TO 2 
BT=" " 

FOR 1=1 TO 15 
GET# 1 ; AT ( I ) 

IF AT ( I ) =CHRT (13) THEN 280 
IF AT(I)=CHRT(10> THEN 340 
IF AT ( I ) =CHRT (32) THEN 340 
NEXT I 

FOR K= 1 TO I 
BT=BT+AT <K) 

NEXT k 
NS=VAL ( BT ) 

PRINT NS 
FOR K=1 TO NS 
BT=" " 

FOR J=1 TO 18 
BT= " " 

FOR 1=1 TO 15 
GET # 1 ; AT ( I ) 

IF AT ( I ) =CHRT (13) THEN 450 
IF AT ( I ) =CHRT ( 32 ) THEN 500 
IF AT ( I ) =CHRT ( 10) THEN 500 
NEXT I 

FOR M= 1 TO I 
BT=BT -t-AT ( M ) 

NEXT M 

DIO (3 ) =VAL (BT) 

DIO ( J) = (DIO(J> -DC (J) ) /1000. 

IF K= 1 THEN REF=DIO (18) 

NEXT J 

FOR J=1 TO 3 
BT= " " 


DIODE SENSITIVITY" 



5 GO 
590 
600 
610 
620 
630 
640 
650 
660 
665 
670 
680 
690 
695 
700 
705 
710 
720 
730 
740 
750 
760 
770 

775 

776 
780 
790 
800 
81 1 
812 
815 
820 
830 
835 

840 

841 

842 

843 

844 

845 
850 
860 
900 
910 
920 
930 
940 
950 
960 
970 
980 
990 

1 000 
1 0 1 0 
1020 
1040 
1 100 
1 1 20 
1 1 30 
1140 
1 150 
1 1 60 
1 200 
1210 
1 220 


GET# 1 : hi ( I > 

IF Af ( I ) =CHR-£ (13) THEN 580 A- 9 

IF A-S ( I ) =CHR-£ ( 32 ) THEN 630 
IF AT < I ) =CHR£ (10) THEN 630 
NEXT I 

FOR M= 1 TO I 
bt=--bt+at (M) 

NEXT M 

ANG (J ) =VAL (B4 ) 

NEXT J 

FOR 1=1 TO 15 

ALF'HA= ANG ( 2 1 - ANG ( 1 ) 

X=SIN (DIANG ( I ) *DTRD) *COS (AUPHA*DTRD) 

NAD= FN ACOS ( X ) 

L=L+ 1 

NADANG (L) =NAD 

COMF- < D I ANG ( I > -90 ) *DTRD 

DEN= ( B I N ( NADANG ( L ) ) *S I N ( ALF'H A*DTRD ) ) 

IF DEN=0. THEN 750 

X= (COS (COMP) -COS (NADANG (L) ) *COS < ALPHA*DTRD) ) /DEN 
IF DEN=0. THEN X=. 00001 
IF DIANG (1X90. THEN X=-X 

IF NADANG (L) T. 000001 THEN 775: ELSE X=1 . 0 
IF X :• 1 . THEN X= 1 . 

IF X- -1. THEN X=- 1 . 

RAZANG (L) = FN ACOS ( X ) 

IF ALPHA ,0. THEN RAZANG (L) =RAZANG (L) +180*DTRD 
IF RAZANG (LKO. THEN RAZANG (L) =360. *DTRD+RAZANG <L) 
MEAS (L , LL) =DIO ( I ) /SENS ( I ) *REF/DIO< 18) 

VOLT ( L ) =ME AS ( L , 2 ) — MEAS ( L , 1 ) 

PRINT L „ NADANG (L) /DTRD, RAZANG (L> /DTRD , VOLT (L) 

NEXT I 
PR I NT: PRINT 
AVZEN=AVZEN+ANG ( 1 ) /NS 
NEXT k 
CLOSE# 1 

IF LL=2 THEN 850 

DPEN#1 AS INPUT , " . D2/FIN60. Z60CUBE" 

L=0: AVZEN=0. 

NEXT LL 
NMEAS=L 
GOSUB 1200 

REM SUBROUTINE INTERP 
DELTA=30*DTRD 
SUMF=0. : SUMW=0. 

FOR 1=1 TO NMEAS 

Z=NADANG < I ) : A=RAZ ANG ( I ) 

IF ABS (Z-ZINT) ,-DELTA THEN 1040 
IF ABS (A-AINT) /-DELTA THEN 1040 

THET A=COS ( Z ) *C03(ZINT) +SIN ( Z ) *SIN ( Z I NT) *COS (A-AINT) 

IF THETA/ 0. THEN 1000: ELSE 1000 
THETA= FN ACOS (THETA) 

GOSUB 1100 
SUMW=SUMW+W 
SUMF=SUMF+W*VOLT ( I ) 

NEXT I 

REM SUBROUTINE WEIGHT 
WGHT= 1 . E+06 

IF THETA '.0001 THEN 1140: ELSE 1150 
W=WGHT: RETURN 
WGHT =DELTA/ THETA 
W=WGHT: RETURN 
REM SUBROUTINE MAX- 
MAX=0 
M T M - 1 on 


MIN 



•1237 

1240 
1 250 
1260 
1265 
1 300 
1 3 1 0 
1315 
1320 
1325 
1 330 
1335 
1340 
1345 
1350 
1355 
1360 
1365 
1370 
1375 
1380 
1390 
1395 
1400 
1410 
1420 
1430 

1435 

1436 

1437 

1438 

1439 

1440 

1441 
1445 
1 500 
1510 
1520 
1530 
1540 
1550 
1560 
1570 
1580 
1590 
1 600 
1610 
1620 
1630 
1640 
1650 
1 660 
1 670 
1680 
1690 
1700 
17 10 
1720 
1730 
1740 
1750 
1760 
1770 

1 7 PC) 


IF N A DANG ( I ) 60*DTRD THEN 1260 

IF VOLT ( 1 1 MAX THEN MAX=VOLT < I ) A _ 10 

IF VOLT (I) MIN THEN MIN=VOLT(I) 

NEXT I 

PRINT MIN, MAX 

REM SUBROUTINE INTEGRATE RADIANCES 
SUM=0. : T 0TNUM=0 

FOR J = 1 TO 36 PMWAL PAGE as 

NR ( J ) =0 : NRN ( J ) =0 QE POOR QUALITY 

NEXT J 

NINC=90*DTRD/36. 

FOR 1=1 TO NMEAS 

J = INT (NADANG ( I ) /NINO+l 
NR ( J ) =NR ( J ) +VOLT < I ) 

NRN ( J > =NRN ( J ) + 1 
NEXT I 

FOR 1=1 TO 36 

IF NRN ( I ) =0 THEN 1395 
ANG= (2*1-1 ) /2*NINC 

SUM=SUM+NR ( I ) /NRN < I ) *COS (ANG) *SIN ( ANG) *NINC 
TOTNUM=TQTNUM+NRN ( I ) 

N I NC=90 . *DTRD/36 . : GOTO 1400 
N I NC=N I NC+N I NC 
NEXT I 

FLUXUF'=3. 14159+2*SUM 

PRINT FLUX UP 

STOP 

0PEN#5 AS OUTPUT, ".SI LENTYPE" 

OUTPUT #5 

FOR F=8 TO 375 STEP 15 

VOLT (K ) = (ME AS <K ,2) -MEAS (K, 1>*( 1-F* ( 1+TAN (NADANG (K> > > > > /F / ( 1+TAN <NADAN 
G (K) ) ) 

PRINT MEAS (K , 1 ) ,MEAS<K,2) ,VOLT(K) , NADANG (K) /DTRD , ( 1 -COS (NADANG (K) ) ) 
NEXT K 
PRINT: PRINT 
END 

0PEN#2 , " . GRAF IX" 

INVOKE" . Dl/BGRAF. INV" 

PERFORM INITGRAFIX 
PERFORM GRAFIXMODE (7.1 ,7.1) 

PERFORM FILLPORT 

RAD=60 

GOSUB 1880 

PERFORM PENCOLOR ( 7. 1 5 ) 

NF'0INTS=72 
PERFORM GRAFIXON 
FOR H =1 TO 3 
X1=1>30*DTRD 
R=RAD*SIN(X1) 

PERFORM MO VETO <7. (R* 140/ 192*2+ 140) ,7.96) 

FOR 1=1 TO 365 STEP 5 

X=R*COS ( I*DTRD) *140/192+70: X=2*X 
Y=Rx-SIN ( I*DTRD) *192/140+96 
PERFORM LINETO (7.X , 7.Y) 

NEXT I 
NEXT K 

ANGI NC=30*DTRD 
PERFORM MOVETO (7.140, 7.96) 

FOR 1=1 TO 12 
ANG= ANG I NC * I 

X=RAD*COS (ANG) *140/ 192+70: X=2*X 
Y=RAD*SIN ( ANG) *192/140+96 
PERFORM LINETO (7.X ,7. Y) 

PERFORM MOVETO (7.140, 7.96) 

MCTY-r T 



A- 1 1 


4 7 r v0 
1995 
1800 
1802 
1804 
1806 
1 808 
1810 
1812 
1814 
1816 
1818 
1820 
1822 
1825 
1 830 
1840 
1850 
1860 
1870 
1875 
1 880 
1885 
1890 
1895 
1 900 
1905 
1910 
1 920 
1925 
1930 
1935 
1940 
1950 
1960 
1970 
1980 
1 990 
2000 
20 1 0 

2020 

er 

.t- 

2030 
2040 
2050 
3000 
30 1 0 
3020 
3030 
3040 
3050 
3060 
3070 
3 1 00 
3110 
3 1 20 
3 1 30 
3140 
3 1 50 
3155 
3160 
3170 
3180 
3185 

3 1 90 


PERFORM MO VETO < 7.230 , 7.98 > 

PR I NT# 2 ; " 1 80 11 
PERFORM MOVETO <7.10,7.188) 

PR I NT# 2 USING 1R04;AVZEN 
IMAGE " Z = “ , ##. 

PERFORM MOVETO (7.180,7.188) 

PRINT# 2 USING 1810;FLUXUP 
IMAGE “F=“,#.## 

PERFORM MO VETO (7.10,7.8) 

PR I NT#2 USING 1816; MIN 
IMAGE “MI N= “ , # . # 

PERFORM MOVETO (7.180,7.8) 

PR I NT #2 USING 1822; MAX 
IMAGE “MAX= U ,#.# 

PERFORM MOVETO (7.126,7.188) 

PR I NT #2 ; “90“ 

PERFORM MOVETO (7.30,7.98) 

PR I NT #2 ; " 0 " 

PERFORM MOVETO (7.120, 7.8) 

PR I NT#2 ; "270" 

GOTO 3000 
PERFORM GRAFIXON 
PERFORM INITGRAFIX 
PERFORM X F PORT I ON ( 7.0 ) 

CQLINC= (MAX -MIN) 710. 

FOR K=1 TO NMEAS 

IF NADANG (K) >60*DTRD THEN 2025 
R=-RAD*S I N ( NADANG ( K ) ) 

X=R*COS ( RAZANG (K) ) *140/ 192+70: X = X*2 
XS=3*C0S (NADANG <K> ) +1 : XS=XS*2 
Y=R*S IN (RAZANG <K) ) *192/140+96 
YS=6*C0S (NADANG (K) ) +3 
IF COLINOO. THEN COLOR* 10: ELSE COLOR* (VOLT <K) — MIN) /COL INC 
COLOR* I NT (COLOR) 

PERFORM F I LLCOLOR ( 7.C0L0R ) 

R*X+XS 

L=X-XS 

T=Y+YS 

b=y-ys 

PERFORM V I EWPORT ( 7.L , 7.R , 7.B , 7.T ) 

PERFORM FILLPORT 
NEXT K 

PERFORM INITGRAFIX 

RETURN 

END 

REM SUBROUTINE PRINCIPAL PLANE RADIANCE GRAPH 

REM FIND THE MAX AND MIN VOLT IN THE PRINCIPAL PLANE 

MAX*0 . ; M I N= 1 00 - 

FOR 1=8 TO (NS-1 ) *15+8 STEP 15 

IF VOLT<I>:;MAX THEN MAX=VOLT(I) 

IF VOLT ( I ) <MIN THEN MIN=VOLT(I) 

NEXT I 

PRINT MAX, MIN 

REM SET XSCALE AND XINC 

XRANGE=NADANG <8) +NADANG < (NS-1 ) *15+8) 

XRANGE=X RANGE /DTRD 

XRANGE= I NT ( ( XRANGE/2+ 10) / 10) 

X I NC=40 

REM SET YSCALE AND YINC 
MID=INT ( (NS+1 ) /2) 

FOR 1=8 TO (NS-1 ) *15+8 STEP 15 

IF NADANG ( I ) <-'8*DTRD THEN NORM=VOLT ( I ) 

NEXT I 
J=0 ; 

pop t=p rnfM^-nns+R rtfp is 



.32 1 0 NR < J > =7GLT < I ) NORM 

3220 NEXT I 

3230 M A X =M AX / NORM : M I N=M IN/ NORM A-12 . 

3240 REM SET Y 'RANGE AND YINC 

3250 YINC =15 

3260 YRANGE=MAX — M I N 

3265 PERFORM GRAF I X MODE ( 7.2 , 7.2 ) 

3270 PERFORM GRAF IX ON 

3275 PERFORM FILLPORT 

3280 PERFORM MOVETO ( 7.280 , 7.30 ) 

3290 PERFORM MOVEREL <7.-3 , 7.3) 

33 1 0 PERFORM MOVETO ( 7.280 , 7.30 ) 

3320 PERFORM MOVEREL < 7.0 , 7.- 1 0 ) 

3330 ANG=0 

3340 PR I NT #2 ANG 

3350 FOR K=- 1 TO 1 STEP 2 

3360 FOR 1=1 TO X RANGE 

3370 XR=280+K * I *40: YR=30 

3380 PERFORM MOVETO <7.XR-1 , 7.YR+3) 

3400 PRINT#2; " I " 

3410 PERFORM MOVETO < 7.XR , 7.YR ) 

3420 PERFORM MOVEREL < 7.-3 , 7.- 10) 

3430 ANG= 1*10 

3440 PRINT#2:ANG 

3450 NEXT I 

3470 NEXT K 

3490 LL=280-XRANGE*40 

3500 LR= 280+ X RANGE *40 

3510 PERFORM MOVETO <7.LL , 7.30 ) 

3520 PERFORM L INETO < V.LR , 7.30 ) 

3530 PERFORM MOVETO ( 7.280 , 7.30 ) 

3540 FOR J = 1 TO 3 

3550 X R = 2 B 0 : YR=30+J*5O 

3560 PERFORM MOVETO (7.XR-3 , 7.YR+2) 

3570 PRINT#2; 

3575 PERFORM MOVETO (7.XR , 7.YR) 

3580 PERFORM MOVEREL ( 7. 1 0 , 7.0 ) 

3590 PRINT#2;J 

3600 NEXT J 

3610 PERFORM MOVETO < 7.XR , 7.YR ) 

3620 PERFORM L INETO < 7.280 , 7.30) 

3622 PERFORM MOVETO < 7.280 , 7.30 ) 

3625 J=0: L=0 

3630 FOR 1=8 TO <NS-1 ) *15+8 STEP 15 
3635 IF VOLT < I ) =0 THEN 3710 

3640 k= 1 

3645 J=J+1 

3650 IF RAZANG(I) .'5*DTRD THEN K=-l 

3660 X=280+F +NADANG < I ) /DTRD*4 

3670 Y=NR ( J ) *50+30 

3675 PRINT NR ( J ) , NADANG < I ) /DTRD , RAZANG < I ) /DTRD , X , Y 

3680 X=INT <X ) 

3690 Y=INT <Y) 

3695 IF ABS < X-509) <3 THEN 3710 

3700 IF L=0 THEN PERFORM MOVETO <7.X ,7.Y) : ELBE PERFORM L INETO (7.X ,7.Y) 

3701 L=L+ 1 

3705 PERFORM DOT AT (7.X , 7.Y) 

3710 NEXT I 

3720 PERFORM MOVETO < 7. 1 00 , 7. 1 80 ) 

3730 PRINT#2 USING 3740; TT 

3740 IMAGE "RELATIVE RADIANCES IN THE PRINCIPAL PLANE FOR " , 1 

OA 

3750 PERFORM MOVETO <7.40 , 7. 170) 

3760 F'R I NT #2 USING 3770;AVZEN 

3770 tm&hp "qm ar 7fnttm amgi F= " .##. . "ofr" 



3790 . 

3800 

3810 

3820 

3830 

3840 

1 1438 


PR I NT #2 " ANT I -SOLAR S I DE " 

PERFORM MOVETO < 7.400 , 7. 1 0 ) 

PR I NT 42 " SOLAR SIDE" A- 13 

PERFORM MOVETO (7.240,7.10) 

PRINT#2" NADIR ANGLE" 

END 

VOLT (K ) = (ME AS (H',2) -MEAS (K, 1 > *< 1-F* ( 1+TAN (NADANG(R') ) ) ) ) / (F* ( 1+TAN (NADAN 
G ( K ) ) ) ) 



A- 14 


10 REM PROGRAM CFOS 

20 PRINT" CFOS CALIBRATION PROGRAM" 

40 DIM DI ANG (50) 

45 DIM DARK' (500) 

50 DIM DIO (40) ,ANG<3) , NR(36) , NRN(36) 

55 DIM AT (15) 

56 DIM VALUE (500,20) 

57 DIM DC (19) 

58 DIM SENS (10, 16) 

60 DATA 20 , 30 , 40 , 50 , 60 , 70 , 80 , 90 , 1 00 ,110,1 20 , 1 30 ,140,150,1 60 

70 DATA 1 . 00 , . 9 1 , 1 . 07 , . 98 , . 99 , . 96 , 1 . 2 , 1 . 09 , 1 . 03 , 1 . 23 , . 89 , 1 . 02 , 1 . 1 0 , . 88 , . 96' 
75 DIM BDR (50) 

190 DTRD=3. 14159/180. 

200 REM DEFINE ARCOS FUNCTION 

210 DEF FN ACOS ( X ) =ATN (SQR ( 1 -X*X ) /X ) 

220 L=0 

225 AVZEN=0 

250 OF'EN# 1 AS I NF'UT ," /AREA1 /CFOS. BACK/D. 07 1 684 . 1244" 

251 TT=" TEST. DARK" 

260 B-£= " " 

280 INPUT# 1 , 1 ; BT 
370 NS=VAL(BT) 

375 NS=NS-2 
380 PRINT NS 
400 FOR Y — 1 TO NS 
410 BT=" " 

420 I NF’UT # 1 , K+ 1 ; BT 

430 PGS1=1 • 

435 FOR J=1 TO 20 

440 SPACE=F'OS 1 

450 F'0S2=INSTR (BT , " ", SPACE) 

460 N0C=F'QS2-P0S1 

470 VALUE (K , J ) =VAL (MIDT ( , F'OSl , NOC) ) 

480 F'OSl =F'0S2+1 

490 PRINT VALUE (K,J) 

500 NEXT J 

510 NEXT K 

520 REF=VALUE (1,17) 

530 NS2=NS/2 
540 FOR J=1 TO NS 

550 VALUE ( J ,8) =VALUE ( J ,8) *REF/ VALUE ( J , 17) 

560 DARK <J)=VALUE(J,9) *0.761652+163.5 

580 NEXT J 

600 F=. 01 

610 FOR J=1 TO NS 

620 DIANG ( J ) = VALUE ( J , 19) 

630 IF DIANG ( J ) <0. THEN D I ANG ( J ) =360. -DI ANG ( J ) 

640 PRINT DIANG (J) 

650 NEXT J 

660 FOR J=1 TO NS 

670 DIANG (J)=( (DIANG ( J ) +VALUE ( J , 18) >-360)*DTRD 

680 DIANG (J)=ABS (DIANG (J) ) 

685 SINE=ABS (SIN (DIANG ( J) ) ) 

686 COSE=ABS(COS (DIANG (J) ) ) 

690 DIO ( J ) = (VALUE ( J ,8) —DARK ( J ) ) 

700 NEXT J 

701 DI0(7)=(DI0(8> +DI0(6) >/2. 

710 FOR J=1 TO NS 

720 PRINT DIO (J ) /DIO (7) , D I ANG ( J ) /DTRD 

730 NEXT J 

740 SUM=0. 



753 SINE=ABS(SIN<DIANGO.) ) ) 

756 COSE—ABS ( COS ( D I ANG < h ) ) ) A_15 

760 SUM=SUM+DIO (H ) *COSE*S I NE* 1 0*DTRD* 1 *3 . 14159 

780 NEXT K 

790 SOL ANG= ( 1 . O-COS ( D I ANG ( 1 ) ) ) *2 . 0*3 .14159 

795 0PEN#2 AS OUTPUT , . S ILENTYPE" 

796 0UTPUT#2 

800 FOR h = 1 TO NS 

810 BDR (F) =DIO <K) *3. 14159/SUM 

820 PRINT BDR <K> ,COS(DIANG <K) ) 

830 NEXT K 

840 END 

850 NMEAS=L 

900 FOR K=1 TO 5 

910 FOR J=1 TO 18 

920 DC ( J ) =DC ( J ) +VALUE (K , J ) /5 

930 NEXT J 

940 NEXT K 

950 FOR K = 1 TO 19 

960 PRINT DC ( k ) 

970 NEXT K 

1000 FOR W = 6 TO NS 

1010 J = (K-5) -15* I NT ( (K--5) / 15) 

1020 IF J=0 THEN J=16 

1030 I = I NT ( ( K-6 ) / 1 5 > + 1 

1040 REF= VALUE (5+1 , 1 ) 

1042 SENS (I , J ) = ( VALUE ( K , J ) -DC ( J ) ) /REF*VALUE < 1+5 , 18) /VALUE (K , 18) 

1050 PRINT I , J , K , VALUE < k , J ) , SENS < I , J ) 

1060 NEXT K 

1100 END 

1200 FOR 1=1 TO 4: FOR J=2 TO 16: SENS ( I , J) =SENS (I , J ) /SENS ( I , 1 ) : NEXT J:NEXT I 

1210 FOR 1 = 1 TO 5: FOR J = 1 TO 16:F'RINT SENS < I , J ) : NEXT J:NEXT I 

1300 0F'EN#3 AS OUTPUT , " . S ILENTYPE" 

1310 OUTPUT #3 

1320 FOR J=1 TO 16 

1330 PRINT SENS ( 1 , J ) ,SENS(2,J) , SENS (3, J) , SENS <4 , J ) ,SENS(5,J> 

1340 NEXT J 


! 
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REM PROGRAM CFOS 

PRINT" CFOS CALIBRATION PROGRAM" 

DIM D I ANG ( 50 ) 

DIM DARK (500) 

DIM DIO (40) , ANG (3) ,NR(36> ,NRN(36) 

DIM AT (15) 

DIM VALUE (500,20) 

DIM DC ( 1 9 ) 

DIM SENS (10, 16) 

DAT A 20 , 30 , 40 , 50 , 60 , 70 , 80 , 90 , 1 00 ,110,120,130,140,150,1 60 
DATA 1 . 00 ,.91,1. 07 , . 98 , . 99 , . 96 , 1 . 2 , 1 . 09 , 1 . 03 , 1 . 23 , . 89 , 1 . 02 
DIM BDR ( 50 ) 

DTRD=3. 14159/180. 

REM DEFINE ARCOS FUNCTION 

DEF FN ACOS ( X ) =ATN (SQR <1-X*X)/X) 

L=0 

AVZEN=0 

OPEN# 1 AS INPUT, "/AREA 1 /CFOS. OFF/D. 07 1984. 0826" 

TT=" TEST. DARK " 

BT=" " 

INF'UT#1 , 1 ; B$ 

NS=VAL (BT> 

NS=NS-2 
PRINT NS 
FOR K= 1 TO NS 
BT=" " 

INPUT# 1 , K+l ; BT 
P0S1=1 

FOR J=1 TO 20 
SPACE=F'0S1 

F’0S2= I NSTR ( BT , " ", SPACE) 

N0C=P0S2-F’0S 1 

VALUE (K , J ) =VAL ( M I D$ ( BT , POS 1 , NOC ) ) 

P0S1=F'0S2+1 
PRINT VALUE <K,J) 

NEXT J 
NEXT K 

REF=VALUE ( 1 , 17) 

NS2=NS/2 
FOR J=1 TO NS 

VALUE ( J , 8) = VALUE ( J , 8) *REF/ VALUE (J ,17) 

DARK ( J ) =VALUE ( J , 9 ) *0 . 87 1 844+47 . 374 
NEXT J 
F = . 0 1 

FOR J=1 TO NS 

D I ANG ( J ) =VALUE ( J , 19) 

IF DIANG(J) (0. THEN DIANG ( J ) =360. -DI ANG ( J ) 

PRINT DIANG (J) 

NEXT J 

FOR J=1 TO NS 

DIANG (J) = ( ( D I ANG < J ) +VALUE < J , 1 8 ) >-360>*DTRD . 

DIANG (J)=ABS (DIANG (J) ) 

S I NE= ABS ( S I N ( D I ANG ( J ) ) ) 

COSE=ABS (COS (DIANG ( J ) ) ) 

DICK J) = < VALUE (J, 8) -DARK (J) ) 

NEXT J 

DIO (7) = (DI0(6)+DI0(8) ) /2. 

FOR J=1 TO NS 

PRINT DIO(J) /DIO (7) , DIANG (J) /DTRD 
NEXT J 
SI = . 


10, . 88, . 96 



755 SINE=ABS <SIN<DIANG<IO ) ) 

756 C0SE=ABS <COS<DIANG<I< ) > ) A17 

760 5UM=SUM+DI0 O.) *COSE*SINE*10*DTRD*1*3. 14159 

73 0 NEXT F 

790 SOL ANG= < 1 . 0— COS < D I ANG < 1 ) ) ) *2.0*3. 14159 

795 0F‘EN#2 AS OUTPUT SILENTYF'E" 

796 0UTPUT#2 

800 FOR H = 1 TO NS 

810 BDR (1 ) =DIO <k) *3. 14159/SUN 

820 PRINT BDR <F) ,COS (DIANG (k) ) 

830 NEXT V 

840 END 

850 NMEAS=L 

900 FOR k=l TO 5 

910 FOR J = 1 TO 18 

920 DC < J ) =DC < J ) + VALUE <K ,J ) /5 

930 NEXT J 

940 NEXT K 

950 FOR K= 1 TO 19 

960 PRINT DC <K) 

970 NEXT K 

1000 FOR K=6 TO NS 

1010 J= <K-5) -15+INT < <K— 5) /15) 

1020 IF J=0 THEN J = 16 

1030 I = I NT < < K— 6 ) / 1 5) + 1 

1040 REF=VALUE<5+I , 1) 

1042 SENS < I , J) = < VALUE <K, J) -DC (J) ) /REF* VALUE <1+5,18) /VALUE <K, 18) 

1050 PRINT I , J , k , VALUE <K , J ) , SENS < I , J ) 

1060 NEXT K 

1 1 00 END 

1200 FOR 1=1 TO 4 : FOR J=2 TO 16: SENS ( I , J ) =SENS < I , J ) /SENS < I , 1 ) : NEXT J:NEXT I 

1210 FOR 1=1 TO 5: FOR J=1 TO 16: PRINT SENS ( I , J ) : NEXT J:NEXT I 

1300 QPEN#3 AS OUTPUT SILENTYPE" 

1310 0UTPUT#3 
1320 FOR J=1 TO 16 

1330 PRINT SENS < 1 , J ) , SENS <2 , J ) , SENS (3 , J ) , SENS <4 , J ) , SENS <5 , J ) 

1340 NEXT J 
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10 REM PROGRAM CFOS 

20 PRINT" CFOS BI-DIRECTIONAL REFLECTANCE PROGRAM 

40 DIM D I ANG < 15) , SENS ( 15) , COFFS (15), NR ( 50) , NRN (50) 

45 DIM DART (20, 15) , DIO (20, 17) ,BDR (20, 15) , OFFS ( 15) ,SLOP (15) 
50 DIM ANG (20, 3) ,NAD(20, 15) ,AZI (20, 15) 

60 DATA 20 , 30 , 40 , 50 , 60 , 70 , 80 , 90 , 1 00 ,110,1 20 , 1 30 , 1 40 , 1 50 , 1 60 
80 FOR 1=1 TO 15 
90 READ DIANG(I) 

100 NEXT I 

110 OPEN# 1 AS INPUT, "/AREA 1 /CFOS. CAL/COEFF" 

120 0PEN#2 AS INPUT, "/AREA1 /CFOS. CAL/ADJUST" 

130 0PEN#3 AS INPUT, "/AREAl /CFOS. CAL/SENSIV" 

140 FOR J = 1 TO 15 

150 READ# 1 , J;SLOP(J) ,OFFS(J) 

160 READ#2, J; COFFS (J) 

170 READ#3 , J : SENS ( J ) 

180 OFFS ( J ) =COFFS ( J ) 

190 NEXT J 

200 CLOSE 

205 DTRD=3. 14159/180. 

210 REM DEFINE ARCOS FUNCTION 

215 DEF FN ACOS ( X ) =ATN ( SQR ( 1— X*X ) /X ) 

220 PRINT" INPUT DATA FILE NAME " 

225 INPUT Df 

250 OPEN# 1 AS INPUT, "/AREA 1 /CFOS. BACK/" +D4? 

260 B*=”" 

280 INF'UT#1 , 1 ; 

370 NS=VAL(BT) 

375 NS=NS— 2 

380 PRINT NS 

400 FOR K = 1 TO NS 

410 B*=" " 

420 INPUTttl ,K+1;B4? 

430 F'OS 1 = 1 

435 FOR J= 3 TO 17 

440 SF’ACE=POS 1 

445 P0S2=INSTR<B$, " ", SPACE) 

450 N0C=P0S2-P0S1 

455 D 1 0 ( K , J ) =OAL ( M I D$ ( B# , F'OS 1 ,NOC) ) 

460 P0S1=P0S2+1 

465 PRINT DIO (K , J) 

470 NEXT J 

475 FOR J= 1 8 TO 20 

480 SF’ACE=POS 1 

485 P0S2=INSTR(B$, " ", SPACE) 

490 N0C=P0S2-P0S1 

4P5 ANG (K , J-l 7 ) =VAL < MID£ (B$ , F'OSl , NOC) ) 

500 P0S1=P0S2+1 

505 PRINT ANG (K, J-l 7) 

510 NEXT J 

515 NEXT Y 

520 REF=DI0( 1 , 17) 

530 NS2=NS/2 

540 FOR J=1 TO NS 

545 FOR K= 1 TO 15 

550 DARK ( J , K ) =D 1 0 ( J , 1 6 ) *SLOF' ( K ) +OFFS < K ) 

560 DIO ( J ,1 ) = (DIO ( J , (■ ) -DART ( J , K ) ) /SENS(K ) 

570 DIO(J ,K)=DIO(J,K)*REF/DIO(J, 17) 

590 NEXT K 

600 NEXT J 

610 FOR J=1 TO NS 


1 ' 


.635 ■ ANG ( J , 2 ) --ANG ( J , 2 > 

640 PRINT ANG (3,2) 

650 FOR 1=1 TO 15 A- 19 ORIGINAL PAGE fs 

680 ALPHA= ANG ( J , 2 ) — ANG ( J , 1 ) OE POOR QUALITY 

690 X = SIN (DI ANG (K ) *DTRD) *COS ( ALF'HA*DTRD) 

700 NAD ( J , k) = FN ACOS ( X ) 

710 CQHP= (DIANG (k) -90) *DTRD 

720 DEN= ( S I N ( NAD ( J , K ) > *SIN < ALPHA*DTRD) ) 

730 IF DEN=0 . THEN 750 

740 X= (COS (COMP) -COS (NAD < J ,K> ) *COS < ALPHA*DTRD) ) /DEN 

750 IF DEN=0. THEN X=0. 00001 

760 IF DI ANG (h ) <'90. THEN X=-X 

770 IF NAD ( J, K ) .0.00001 THEN 775:ELSE X=1.0 

775 IF X ' 1 . THEN X=1 . 

776 IF X / — 1 . THEN X=-1.0 

780 AZ I < J , 1 ) = FN ACOS ( X ) 

790 I F ALPHA/ 0 . 0 THEN AZ I < J , K ) =AZ I < J , K ) + 1 80 . 0*DTRD 

800 IF AZI(J,KK0. THEN AZ I ( J , K > =360 . 0*DTRD+AZ I ( J , K) 

810 PRINT J , K , NAD ( J , K) /DTRD , AZ I ( J , K) /DTRD , DIO ( J ,K) 

820 NEXT k 

830 NEXT J 

840 SUM=0 

850 T0TNUM=0 

860 FOR J=1 TO 36 

870 NR ( J ) =0 

880 NRN ( J ) =0 

900 NEXT J 

910 NINC=90*DTRD/36. 

920 FOR J = 1 TO 13 
930 FDR 1=1 TO 15 

940 i=inj (NAD (J,K> /NINC) +1 

950 NR ( I ) =NR ( I ) +D I 0 ( J , K ) 

960 NRN ( I ) =NRN ( I ) + 1 

970 NEXT k 

980 NEXT J 

1000 FOR 1=1 TO 36 

1010 IF NRN ( I ) =0 THEN 1070 

1020 ANG= (2*1-1 > /2*NINC 

1030 SUM=SUM+NR ( I ) /NRN ( I ) *COS ( ANG) *SIN < ANG) *NINC 

1 040 TOTNUM=TQTNUM+NRN ( I ) 

1050 NINC=90. *DTRD/36. 

1060 GOTO 1090 

1070 N I NC=N I NC+N INC 

1090 NEXT I 

1100 FLUXUP=3. 14159*SUM*2. 

1110 PR I NT " UPWARD FLUX DENSITY" ; FLUXUP 

1120 NAX=0 

1130 NIN=90 

1140 FOR J=1 TO 13 

1150 FOR K=1 TO 15 

1160 IF NAD ( J , K) ./-MAX THEN MAX=NAD(J,K) 

1170 IF NAD<J ,KXMIN THEN MIN=NAD(J,K> 

1180 NEXT K 

1190 NEXT J 

1200 PRINT N IN, MAX 

1210 S0LANG=2.0*3. 14159* ( 1 . -COS (MAX ) ) 

1220 FRACSA=SOLANG/ (2. *3. 14159) 

1230 PRINT SOLANG , FRACSA 
1300 FOR J=1 TO 13 
1310 FOR y = 1 TO 15 

1320 BDR ( J , Hi ) =D 1 0 < J , K ) *3 . 14159/FLUXUP 

1330 PRINT J , k , NAD < J , K) /DTRD , AZ I < J , K) /DTRD , BDR ( J , K) 

1340 NEXT K 

1350 NEXT J 

1 "’An MAYr 
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1375 
1380 
1390 
1400 
1410 
1480 
1430 
1 440 
1460 
1470 
1475 
1480 
1490 
1 500 
1510 
1520 
1 530 
1540 
1 550 


1 = 1 

FOR J=1 TO 13 
FOR 1=1 TO 15 

IF BDR ( J ,10 /MAX THEN MAX=BDR < J ,1- ) 

IF BDR (J ,F ) - MIN THEN MIN=BDR<J,K) 

NEXT K 
NEXT J 

PRINT MIN, MAX 

0PEN#3 AS OUTPUT, "/AREA1 /CFOS. BDR/"+D* 

WR I TE#3 , 1 ; NS , FLUX UP , M IN , MAX 
1 = 1 

FOR J=1 TO NS 
FOR F=1 TO 15 
1 = 1 + 1 

WRITE#3 , I ; J , K , NAD ( J , R) , A2I (J,K) ,BDR(J ,K) 
NEXT K 
NEXT J 
CLOSE 
END 
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:i) , XNUM < 2 1 ) , YNUM < 2 1 ) 

, BDA (2,12) , BD ( 30 ) , N A ( 1 4 ) 


234 

235 

236 

237 

239 

240 
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247 

248 

249 

250 

251 

'-icr 

v_> 

253 

254 

255 

256 


HOME 

X0=0 . 0 : Y0=0 . 0 : X 1 =6 : Y 1 =4 
DIM V(ll) ,SCALEX(21) ,SCALEYC 
DIM XG (200) , YG (200> ,BDN (6, 12) 

FOR 1=1 TO 10 

F‘ 1 X = 1 500 : P 1 Y= 1 500 
X 0=0 . 0 : Y0=0 . 0 : X 1 =6 : Y 1 =4 
XMIN=0sYMIN=0. : XMAX= 1 80 : YMAX=3. 0 
NDIVX=12: NDIVY=5: NP=14: DIRUN=1 . 0: DIRISE=0. 

DTRD=3. 14159/180 

SAW=. 1*2. 54: SAH=. 13*2.54 

GOSUB 500 

END 

REM INITIALIZE PLOTER 
0PEN#9 , " . RS232" 

PRINT #9 ; CHR 1 (27) : " . 132; ; 17:": PR I NT #9 ; CHR* ( 27 ) ; " . N; 19; “ ; ’’ IN; " 

F'l X*=CONV* (P IX): P2X*=C0NV* (P2X ) : PI Y*=CONV* (PI Y) : P2Y*=C0NV* (P2Y) 
P2X=P1X+(X1-X0) *1010: P2Y=P1Y+ (Yl-YO) *1010 
IF XO=X 1 THEN P 1 X=500 : P2X= 1 0000 : P 1 Y=600 : P2Y=7500 

P 1 X* =CQNV* ( P 1 X ) :F'2X*=C0NV*(F'2X) : P1Y*=CQNV* (PI Y) : F'2Y*=C0NV* (P2Y) 
PF I ELD*=F‘ 1 X*+ " , "+P1Y*+" , "+P2X*+" , "+P2Y* 

PRINT49; " IP" ; PF IELD*; " 

REM SUBROUTINE TO SCALE AXISES 

REM XO , X 1 , YO , Y 1 DEFINE STARTING AND ENDING POINTS IN INCHES 
REM XM IN , XMAX , YMIN , YMAX ARE THE MAX AND MIN VALUE ON AXISES 
REM SCALEX , SCALEY=SCALED DATA FOR TIC MARKS 
REM NDIVX , NDIVY=NUMBER OF DIVISIONS ON AXISES 
DELTA X= ( X 1 -XO) /NDI VX 
DELT AY* < Yl-YO) /NDI VY 
SCALEX ( 1 ) =X0: SCALEY ( 1 ) =Y0 

FOR I SUB= 1 TO NDI VY: SCALEY < ISUB+1 ) =SCALEY < ISUB) +DELTAY : NEXT 
FOR ISUB=1 TO NDIVX : SCALEX (ISUB+1 )=SCALEX (ISUB) +DELTAX: NEXT 
SC*=CONV* ( XO ) + " , •• +CONV* < X 1 ) + 11 , " +CONV* ( YO ) + " , " +CONV* < Y 1 ) 

PR I NT #9 ; " SC " ; SC* ; " ; " : RETURN 
REM SUBROUTINE TO SET SYMBOL MODE 
SYM 1 *=CONV* ( W ) + " , "+CONV*(H) 

PRINT#9; "SI " ; SYM1*; "SM";SYM*; " ; " 

RETURN 

REM THIS SUBROUTINE PLOTS SYMBOL SYM* AT LOCATION X,Y 
REM PD*0 .FOR PEN UP; PD=1 FOR PEN DOWN; LINET=LINE TYPE(0-4) 

REM LINEL= LENGTH (PERCENT OF DIST. BETWEEN PI AND P2> 

IF PD=0 THEN PRINT#9; "F'U; 11 ; : ELSE F'RINT#9; "F'D; " ; 

IF PD=0 GOTO 239 
IF LINEL=0 THEN LINEL=4 
LT*=CONV* ( L I NET ) + " , " +CONV* ( L I NEL ) 

F'RINT#9; "LT";LT*; " 

SYMB*=SYM* 


SYM2*=C0NV* (X) , "+CONV* (Y) 

F'RINT#9; "PA" ; SYM2* ; " ; " 
RETURN 

REM PLOT AXIS 
PRINT49; "F'U; " 

X0*=C0NV* (XO) : YO*=CONV* ( YO) 
PRINT#9; "PA" ; XO*; " , " ; YO*; " 
PRINT49; "PA" ; XO*; " , " ; YO*; " 
PR I NT #9 ; " TL 1 . 0 , 0 . O ; PA " ; XO* 
FOR I SUB= 1 TO 
PRINT#*; "PA 


X 1 *=CONV* ( X 1 ) : Y 1 *=CONV* ( Y 1 ) 


! "PD; PA" ; X 1 *; " , " ; YO*; " ; " ; "F'U; 
I "PD; PA" ; XO*; " , ";Y1*; "F'U; 

, "YO*; " ; " 

NDIVX+1 : SCALE*=CONV* (SCALEX (ISUB) ) 

; SCALE*; YO*; "; XT; ": FOR 1=1 TO 100: NEXT 
F'R I NT #9 ; " TL 1 . 0 , 0 - 0 ; F’A " ; XO* ; " , " YO* ; " ; " 

FOR I SUB= 1 TO NDIVY+1 : SCALE*=CONV* (SCALEY ( ISUB) ) 

PRINT#9; "PA“; XO* ; SCALE*; ";YT; ": FOR 1=1 TO 100: NEXT 
RETURN 


I: NEXT 


I: NEXT 


ISUB 

ISUB 



A-22 


26 1 , REM LOCATION IS THE LOWER LEFT HAND CORNER 

262 REM DIRUN, DIRISE: SET DIRECTION . 1,0 FOR HORIZONTAL 

263 REM SAW, S AH SET WIDTH AND HEIGHT IN CM. 

264 IF SAW=0 THEN SAW=0. 1 : SAH=0 . 1 5 

266 AXX*=CONV* (AXX) : AYY*=CONV* ( AYY) : AXNUM*=CONV* < AXNUM) 

267 F'R I NT #9 ; "SI" ; CONV* (SAW) ; “ , " ; CONV* (SAH) ; " ; " 

263 PR I NT #9 ; "PA" ; AXX*; " , " ; AYY*; " ; LB" ; AXNUM-?; CHR* (3) 

26R RETURN 

270 REM SUBROUTINE TO PLOT TEXT TXT* AT LOCATION TX X , TYY 

271 REM LOCATION IS THE LOWER LEFT HAND CORNER 

272 REM DIRUN, DIRISE; SET DIRECTION . 1,0 FOR HORIZONTAL 

273 REM SAW , SAH SET WIDTH AND HEIGHT IN CM. 

274 IF SAW=0 THEN SAW=0. 1 : SAH=0. 15 

276 TXX*=CONV* (TXX) : TYY*=CONV* ( TYY ) 

277 PRINT#?; "SI" ; CONV* (SAW) ; " , " ; CONV* (SAH) ; " 

278 PRINT#9; "DI " ;CONV* (DIRUN) CONV* (DIRISE) FOR M=1 TO 999; NEXT M 

279 PRINT#?; "PA" ; TXX*; " , " ; TYY*; " ; LB" ; TXT*;CHR* (3) 

280 RETURN 

281 REM SUBROUTINE TO CENTER DIGITS ON X AXIS 

282 NUMDIGIT=LEN <CQNV£ (AXNUM) ) 

283 AXX=AXX-(NUMDIGIT*SAW) /2/2. 54 

284 RETURN 

285 REM SUBROUTINE TO CENTER DIGITS ON Y AXIS 

236 AYY=AYY-SAH/2/2. 54: RETURN 

287 REM PLOT A CURVE OF NP POINTS IN XG , YG 

283 XF'=(XG(1>— XNUM(5) ) *DX : YP= ( YG < 1 ) -YNUM ( 1 ) ) *DY 

239 SYM2*=C0NV* <XP) +" , "+CONV* ( YP) 

290 PRINT#9; "F’U; F’A" ; SYM2* ; " ; F'D; " 

291 FOR IF'LQT=2 TO NP 

292 XF'= ( XG ( I PLOT) -XNUM ( 5) ) *DX : FOR D=1 TO 200: NEXT D 

293 YF'= ( YG ( I PLOT ) -YNUM ( 1 > ) *D Y 

294 SYM2*=C0NV* ( XP> +" , "+CONV* < YP) 

295 PRINT#?; "PA";SYM2*; " ; " : NEXT I PLOT 

296 F'RINT#9; "F'U; " 

297 RETURN 

299 RETURN 

300 REM SUBROUTINE TO GET ERRORS ON HP 

301 F‘R I NT #9 ; CHR* ( 27 ) ; " .E; " 

302 PRINT#?; "OE; INPUT#?; ERROR*: PRINT ERROR* 

303 RETURN 

500 GOSUB 200 

510 GOSUB 245 

520 FOR 1=1 TO 11 

525 XNUM(I)=<I-1)*0. 10 

530 NEXT I 

535 FOR 1=2 TO 10 

540 XNUM (I+10)=l. 1 - 1 * . 1 

545 NEXT I 

550 FOR 1=1 TO 6 

560 YNUM ( I ) = ( I— 1 ) *. 50 

570 , NEXT I 

575 SAW= . 10*2. 54: SAH=. 10*2.54 

580 AXX=SCALEX (2) 

590 AYY=-. 15*2.54 

600 FOR 1=1 TO NDIVX-1 

610 AXNUM=XNUM (1+5) 

620 GOSUB 281 

625 FOR J=1 TO 100: NEXT J 

630 GOSUB 260 

640 AX X=SCALEX < I +2) 

650 NEXT I 

660 DX= ( X 1-XO) / (XNUM (NDIVX+1 ) —XNUM ( 1 ) ) 

665 DXH=DX 

680 TXT*=" COSINE OF THE OBSERVATION ZENITH" 

6P0 NUMI FT=I FN(TXT*> 
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1150 
1160 
1170 
1 1 80 
1 1 90 
1195 
1200 
1210 
1220 
1225 
1230 
1240 
1250 
1260 


TV V =— . 3*2. 54 
GOSUB 270 

DY=<Y1-Y0) / < YNUM <NDI VY+1 ) — YNUM ( 1 ) ) 
AXX=-. 2*2. 54: AYY=0 
FOR 1=1 TO ND I VY+1 
AXNUM=YNUM(I ) 

GOSUB 285 

FOR J=1 TO 200: NEXT J 
GOSUB 260 

FOR J=1 TO 200: NEXT J 
AYY=SCALEY ( 1+1 ) 

NEXT I 

TXT -£= " NORMAL I ZED RADIANCE" 

TXX=— 0. 35*2.54 
D I RUN=0 . 

DIRISE=1 

NUMLET =LEN ( TXTT ) 

TYY=Yl/2. —NUMLET *SAW/2/2 . 54 
FOR J = 1 TO 800: NEXT J 
GOSUB 270 
FOR 1=1 TO NDIVX+1 
XNUM ( I ) = ( I - 1 ) *0 . 1 
NEXT I 

FOR 1=1 TO NP 
PRINT K 


A-23 


ORIGINAL PAGE sS 
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INPUT BD <K) 

NEXT K 

FOR K= 1 TO NP 

REM PR I NT K : INPUT X G < K ) 

YG<K)=BD<K) 

PRINT X G < K ) , Y G ( K ) 

NEXT K 
DX=5 . 

GOSUB 287 
SYMT=" X " 

W= . 2 : H= . 2 
FOR K= 1 TO NP 

X= ( XG (k) -XNUM (5) ) *DX 
Y= ( YG (K ) -YNUM ( 1 ) ) *DY 
FOR J=1 TO 400: NEXT J 
GOSUB 225 
GOSUB 230 
NEXT K 
END 

TXTT="OCEAN BDR IN THE PRINCIPAL PLANE" 

D I RUN= 1 . 0 : D I R I SE=0 : SYMT= " " : GOSUB 225 

TXX=1 . 25 

TYY=5 

GOSUB 270 

END 

TX X= 1 . 5 
TYY=4 . 5 

SAW=. 075*2. 54: SAH=0. 075*2. 54 

TXTT="0 - COSINE OF THE SOLAR ZENITH = 0." 

GOSUB 270 

END 

TYY=4 . 25 

TXT T= " X - COSINE OF THE SOLAR ZENITH = 0.5" 
GOSUB 270 
END 

TYY=4 . 0 
TXT4?="0 - 0.2 
GOSUB 270 
END 


CLOUD COVER 
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1 

5 

10 
15 
SO 
90 
95 
100 
1 10 
1 15 
120 
125 
150 
200 
201 
202 

203 

204 

205 

206 


HOME 

X 0=0 . O : Y 0=0 . 0 : X 1 =6 : Y 1 =4 

DIM V ( 1 1 ) , SCALE X (21 ) , SCALEY (21) , XNUM(21 ) ,YNUM (21 ) 

DIM XG ( 200 > , YG (200) ,BDN<o, 12) ,BDA<2, 12) ,BD<30> ,NA(14) 

FOR 1=1 TO 10 

p 1 X = 1 500 : P 1 Y= 1 500 
X 0=0 . 0 : Y 0=0 . 0 : X 1 =6 : Y 1 =4 
XM I N=0 : YM I N=0 . : XMAX= 1 SO : YMAX=3 . 0 
NDIVX=12: NYDIV=5: NF-24: DIRUN=1 . 0:DIRISE=0. 

DTRD=3. 14159/180 

SAW= .1*2.54: SAH= . 1 3*2 . 54 

GOSUB 500 

END 

REM INITIALIZE F'LOTER 
OPEN#R, " . RS232" 

F’R I NT #9 ; CHR* ( 27 ) ; " . 132; ; 17: " : F'RINT#9; CHR* (27) ; " . N; 19; " ; " IN; " 

R 1 X *=CONV* ( P 1 X ) :P2X*=C0NV* (P2X) s F'1Y*=C0NV* (F'lY) : F'2Y*=C0NV* (P2Y) : 
P2X=P1X+ (Xl-XO) *1010: P2Y=F'1Y+ (Yl-YO) *1010 
IF XO=X 1 THEN F'1X=500: F'2X=10000: P1Y=600: P2Y=7500 

P1X*=C0NV* (F'l X) :F'2X*=C0NV* (F'2X) : F'l Y*=CONV* (F'lY) s F'2Y*=C0NV* (P2Y) 


207 F'F I ELD*=P 1 X*+" , "+F'1Y*+" , "+F'2X*+" , "+F'2Y* 

208 PRINT#9; " IF'" : PFIELD*; " ; 11 

210 REM SUBROUTINE TO SCALE AXISES 

211 REM X 0 , X 1 , Y 0 , Y 1 DEFINE STARTING AND ENDING POINTS IN INCHES 

212 REM XM I N . XMAX , YM I N , YMAX ARE THE MAX AND MIN VALUE ON AXISES 

213 REM SCALEX , SCALEY =SCALEE> DATA FOR TIC MARKS 

214 REM NDIVX ,NDIVY=NUMBER OF DIVISIONS ON AXISES 

215 DELTAX=(X1-X0) /NDIVX 

216 DELTA Y= (Yl-YO) /NDIVY 


217 

218 

219 

220 
X-X_ 1 



234 

235 

236 

237 

239 

240 

241 

245 

246 

247 

248 

249 

250 

251 


254 


255 

256 


SCALEX .( 1 ) =X0: SCALEY ( 1 ) =Y0 

FOR I SUB= 1 TO NDIVY : SCALEY ( ISUB+1 ) =SCALEY < I SUB ) +DELTAY : NEXT 
FOR I SUB= 1 TO NDIVX : SCALEX (ISUB+1 ) =SCALEX ( ISUB) +DELTAX : NEXT 
SC*=CQNV* (XO) + " , " +CONV* ( X 1 ) + " , "+CONV* (YO) +" , "+CONV* (Y1 ) 

F’R I NT #9 ; " SC " ; SC* ; " ; " : RETURN 

REM SUBROUTINE TO SET SYMBOL MODE 

SYM1*=CQNV* (W) + " , "+CGNV*(H> 

F'R I NT 49 ; "SI" ; SYM1*; "SM" ; SYM$; " ; " 

RETURN 

REM THIS SUBROUTINE PLOTS SYMBOL SYM* AT LOCATION X,Y 

REM F'D=0 FOR PEN UP; PD=1 FOR PEN DOWN; LINET=LINE TYPE (0-4) 

REM LINEL= LENGTH (PERCENT OF DIST. BETWEEN PI AND P2> 

IF PD=0 THEN PRINT #9; "F'U; " ; : ELSE PRINT#9; "PD; " ; 

IF PD=0 GOTO 239 

IF L I NEL=0 THEN LINEL=4 

LT*=CQNV* ( L I NET ) + " , " +CONV* ( L I NEL ) 

F'R I NT #9 ; "LT" ; LT*; "; " 

SYM2*=C0NV* ( X ) + " , " +CONV* (Y) : SYMB*=SYM* 

PRINT49; "F'A";SYM2*; "; " 

RETURN 


REM PLOT AXIS 
F'R I NT #9 ; "F'U; " 


XO*=CONV* (XO) : YO*=CONV* (YO) : X 1 *=CONV* ( X 1 > : Y1*=C0NV* ( Y1 > 


PRINT49; "PA"; XO*; " , " ;Y0*; " 
PRINT#9; "PA" ; XO*; " , " ; YO* ; " 
F'R I NT #9 ; "TL1 . 0,0. 0; F'A" ; XO* 


" ; "F'D; F'A" ; XI*; " , "; YO*; " ; "F'U; " 

" ; "F'D; F'A" ; XO*; " , ";Y1*; " ; " ; "PU; " 

" , "YO*; " ; " 

FOR ISUB=1 TO NDIVX+1 : SCALE*=CONV* (SCALEX ( ISUB) ) 

F'RINTtt9; "F'A" ;SCALE*; " , " ; YO*; " ; XT; " : FOR 1 = 1 TO 100:NEXT I:NEXT 
PRINT49; "TL1 .0,0.0 ; PA" ; XO*; " , "YO*; " ; " 

FOR I SUB= 1 TO ND I V Y+ 1 : SCALE*=CONV* (SCALEY ( ISUB) ) 

PRINT49; "PA" ; XO*; SCALE*; "; YT; ": FOR 1=1 TO 100:NEXT I : NEXT 
RETURN 


ISUB 


ISUB 
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,261 , REM LOCATION IS THE LOWER LEFT HAND CORNER 

062 REM DIRLIN, DIRISE; SET DIRECTION . 1,0 FOR HORIZONTAL 

263 REM SAW , SAH SET WIDTH AND HEIGHT IN CM. 

264 IF SAW=0 THEN SAW=0. 1 ; SAH=0. 15 

266 AXXf=CQNV* (AXX) : AYY*=CON'V* ( AYY ) ; AXNUM*=CONV* ( AXNUM) 

267 PRINT# 1 ?; "SI " ; CONV* (SAW) ; " , " ; CONV* (SAH) ; " ; " 

268 PRINT #9 ; "PA"; AXX*: ", " ; AYY*; " ; LB" ; AXNUM-?; CHR* (3) 

269 RETURN 

270 REM SUBROUTINE TO PLOT TEXT TXT* AT LOCATION TX X , TYY 

271 REM LOCATION IS THE LOWER LEFT HAND CORNER 

272 REM DIRUN, DIRISE; SET DIRECTION . 1,0 FOR HORIZONTAL 

273 REM SAW, SAH SET WIDTH AND HEIGHT IN CM. 

274 IF SAW=0 THEN SAW=0. 1 ; SAH=0. 15 

276 TXX*=CONV* (TXX) : TYY*=CONV* (TYY) 

277 F’RINT#9; "SI" ; CONV* (5AW) ; " , " ; CONV* (SAH) ; " ; " 

278 PR I NT #9 ; " D I " ; CONV* ( D I RUN CON V* ( D I R I SE FOR M=1 TO 999:NEXT M 

279 PRINT#9; "PA"; TXX*; " , " ; TYY*; " ; LB" ; TXT*; CHR* (3) 

280 RETURN 

281 REM SUBROUTINE TO CENTER DIGITS ON X AXIS 

282 NUMD I G I T =LEN ( CONV* (AXNUM) ) 

283 AXX=AXX-(NUMDIGIT*5AW) /2/2.S4 

234 RETURN 

285 REM SUBROUTINE TO CENTER DIGITS ON Y AXIS 

286 AYY=AYY-SAH/2/2. 54: RETURN 

287 REM PLOT A CURVE OF NP POINTS IN XG , YG 

283 XF'= ( XG ( 1 ) -XNUM (5) ) *DX : YF’= ( YG ( 1 ) -YNUM ( 1 ) ) *DY 

289 SYM2*=C0NV*(XP)+" , "+CONV*(YP) 

290 PR I NT#9 ; " F'U ; F'A " ; SYM2* ; " ; PD ; " 

291 FOR IF'LQT=2 TO NF' 

292 XF'= ( XG ( I F’LOT ) -XNUM ( 5) ) *DX : FOR D=1 TO 200:NEXT D 

293 YF'= (YG ( I PLOT) -YNUM ( 1 ) ) *DY 

294 SYM2*=C0NV* ( XP> + " , "+CONV* <YF') 

295 PRINT#9; "PA" ; SYM2*; " ; " : NEXT I PLOT 

296 F'RINT#9; "F'U; " 

297 RETURN 

299 RETURN 

300 REM SUBROUTINE TO GET ERRORS ON HP 

30 1 F'R I NT#9 ; CHR* ( 27 ) ; " . E ; " 

302 PRINT#9; "OE; ": INF'UT#9; ERROR*: PRINT ERROR* 

303 RETURN 

500 GOSUB 200 

510 GOSUB 245 

520 FOR 1=1 TO 11 

525 XNUM < I ) = ( I — 1 ) *0 . 1 0 

530 NEXT I 

535 FOR 1=2 TO 10 

540 XNUM ( I+10)=l . 1-1*. 1 

545 NEXT I 

550 FOR 1=1 TO 6 

560 YNUM ( I ) = ( I - 1 ) * . 50 

570 NEXT I 

575 SAW= . 10*2. 54: SAH=. 10*2.54 

580 AXX=SCALEX (2) 

590 A YY=-. 15*2.54 

600 FOR 1=1 TO NDIVX-1 

610 AXNUM=XNUM( 1+5) 

620 GOSUB 281 

625 FOR J=1 TO 100: NEXT J 

630 GOSUB 260 

640 AXX=SCALEX ( 1+2) 

650 NEXT I 

660 DX= (Xl-XO) / (XNUM(NDIVX+1 ) -XNUM (1 ) ) 

680 TXT*=" COSINE OF THE OBSERVATION ZENITH" 

690 NUMLET=LEN (TXT*) 

700 TX X = X 1 / 2 . -NUMLET + SAW/2 /2 . 54-0 . 30 



,720 , 

GO SUB 270 


-v'-.cr 

DY= (Yl-YO) / (YNUMlNDIVY+l ) -YNUM < 1 ) ) 


726 

AXX=-. 2*2. 54: AYY=0 

• A- 26 

730 

FOR 1=1 TO NDIVY+1 


740 

AXNUN=YNUN ( I ) 


750 

GOSUB 285 


760 

FOR J=1 TO 200: NEXT J 


770 

GOSUB 260 


775 

FOR J=1 TO 200: NEXT J 


780 

AYY=SCALEY ( 1 + 1 ) 


790 

NEXT I 


800 

TX T-f = "NORMAL I ZED RADIANCE" 


810 

TXX=-0. 35*2. 54 


820 

D I RUN=0 . 


830 

D I R I SE= 1 


840 

NUMLET=LEN ( T X TS ) 


850 

TYY=Y 1/2. — NUMLET *SAW/2/2 . 54 


855 

FOR J=1 TO 800: NEXT J 


860 

GOSUB 270 


870 

FOR 1=1 TO NDIVX+1 


880 

X NUN ( I ) = < I - 1 ) *0 . 1 


890 

NEXT I 


900 

FOR K= 1 TO NR 


910 

PRINT K 


920 

REN INPUT BD (K> 


930 

NEXT K 


940 

FOR V=1 TO NP 


950 

XG (K) =1 . +COS < (-57. 5+5*K> *DTRD) 


960 

YG(K)=BD(K) 


970 

PRINT X G ( K ) , Y G < K ) 


980 

NEXT K 


985 

DX=DX*NDIVX/NP 


990 

GOSUB 287 


1000 

SYNT="X" 


1005 

W= . 2 : H= . 2 


1010 

FOR K= 1 TO NP 


1020 

X= ( XG <h ) -XNUM (5) ) *DX 


1 030 

Y= ( YG (K) -YNUN ( 1 ) ) *DY 


1040 

FOR J=1 TO 400: NEXT J 


1045 

GOSUB 225 


1050 

GOSUB 230 


1060 

NEXT K 


1070 

END 


1 100 

TXTT="CFOS BDR IN THE PRINCIPAL PLANE" 


1 1 05 

D I RUN= 1 . 0 : D I R I SE=0 : SYN$= " " : GOSUB 225 


1 1 20 

TX X = 1 . 25 


1 1 30 

TYY=5 


1140 

GOSUB 270 


1145 

END 


1 1 50 

TX X = 1 . 5 


1 1 60 

■ TYY=4 . 5 


1170 

SAW= . 075*2 . 54 : SAH=0 . 075*2 . 54 


1 1 80 

TXT T=" COSINE OF SOLAR ZEN = 1.0" 


1 1 90 

GOSUB 270 


1195 

END 


1200 

TYY=4 . 25 


1210 

TXT " X - CLEAR SURFACE" 


1 220 

GOSUB 270 


1230 

TYY=4. 0 


1240 

TXTT="0 - 0.2 CLOUD COVER" 


1250 

END 
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10 REM PROGRAM CFOS 

20 PRINT" CFOS BI-DIRECTIONAL REFLECTANCE PLOT PROGRAM" 

45 DIM RAD < 1 0 ,10) , RDN ( 30 , 30 > 

50 DIM NXF' < 10 , 10) , NYF' (10,10) , NXD (30) , NYD ( 30) 

60 DATA 2.5,7.5,12.5,17.5,22.5,27.5,32.5,37.5,42.5,47.5,52.5,57.5 

70 DATA 15,45,75,105,135,165,195,225,255,285,315,345 

75 PI =3. 1 4 1 59 : F'0VR2=P 1/2 

80 PR I NT "ENTER THE NUMBER OF X' DATA POINTS AVAILABLE" 

85 INPUT NPX 

90 PR I NT "ENTER THE NUMBER OF Y DATA POINTS AVAILABLE" 

95 INPUT NF'Y 

100 FOR J=1 TO 30 
110 FOR 1=1 TO 30 

120 RDN ( J , K > =0 . 

150 NEXT K 

160 NEXT J 

165 PRINT "ENTER THE NUMBER OF X INTERPOLTION POINTS" 

167 INPUT NX 
170 FOR J=1 TO NX-1 
190 NXD ( J ) =60/NX *J 

200 NEXT J 

216 F'R I NT " ENTER THE NUMBER OF Y INTERPOLATION POINTS" 

217 INPUT NY 

220 FOR J=1 TO NY-1 
230 NYD (J) =60/ NY* J 

240 NEXT J 

250 PR I NT "ENTER THE RADIANCE DATA AS X POS Y POS RAD < X , Y) " 

260 FOR J=1 TO NPX 

270 FOR K= 1 TO NPY 

280 INPUT NXP(J,K) ,NYF'(J,K) ,RAD(J,k) 

290 NEXT K 

300 NEXT J 

400 GOSUB 1000 

500 REM SUBROUTINE TO CALCULATE WEIGHT 

510 DELT A=60/ 1 2 

520 IF THETA ;DELTA GOTO 570 

530 WGHT=1 . 0E+06 

540 IF THETA =.'0.001 THEN RETURN 

550 WGHT =DELTA/ THETA 

560 RETURN 

570 WGHT =0 . 

580 IF (THETA-DELTA) ";=DELTA THEN RETURN 
590 WGHT=2. O-THETA/ DELTA 

600 -RETURN 

700 REM SUBROUTINE TO INTERPOLATE IN BDR ARRAY 
710 SUMF=0. 

720 SUMW=0. 

730 FOR M= 1 TO NPX 

740 FOR N=1 TO NPY 

750 XI=NXP(M,N) 

760 YI=NYP (M , N) 

770 THET A=SOR ( ( X I-XF'RM) * 2 + < Y I-YPRM ) '"2) 

800 ' GOSUB 500 

810 W=WGHT 

820 SUMW=SUMW+W 

830 SUMF=SUMF+W*RAD(M,N) 

840 NEXT N 

850 NEXT M 

860 INTERP=0. 

870 IF SUMW =0. THEN RETURN 

PRO TNTFRP^SI IMF/8UMW 



1 000 
1010 
1020 
1 030 
1040 
1320 
1 330 
1340 
1350 
1 360 
1400 
1410 
1420 
1430 
1440 
1450 
1460 
1470 
1480 
1485 
1490 
1495 
1 500 
1510 
1520 
1 530 
1 535 
1540 
1550 
1560 
1570 


REN CALCULATE RAD VALUES AT REGULAR' ARRAY POINTS 
FOR J=l TO NX-1 
FOR T = 1 TO NY-1 
XPRM=NXD<J1 
YF'RM=NYD < K ) 

GOSUB 700 

RDN ( J , K) =INTERP 

PRINT RDN ( J , L ) 

NEXT K 
NEXT J 

0PEN#4 AS OUTPUT, " / AREA 1 /CFOS. PLOT/" +"STYRO" 

MIN=99999 

MAX=0 

FOR J=1 TO NX-1 
FOR K= 1 TO NY-1 

IF RDN (J,KXMIN THEN MIN=RDN(J,K) 

IF RDN ( J ,K) I; MAX THEN MAX=RDN ( J , K) 

NEXT K 
NEXT J 

PRINT "MIN AND MAX OF INTERPOLATED FIELD ", MIN, MAX 
WRITE#4 , 1 ; NX-1 , NY- 1 , M IN , MAX 
1 = 1 

FOR J=1 TO NX-1 
FOR K= 1 TO NY-1 
1 = 1 + 1 

WR I TE#4 , I ; NXD < J ) , NYD <K) ,RDN<J ,K) 

PRINT NXD(J) ,NYD(K) ,RDN(J,K) 

NEXT K 
NEXT J 
CLOSE 
END 
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20 
40 
50 
60 
70 
80 
1 00 
1 10 
115 
120 
130 
200 
201 
202 

203 

204 

205 

206 

207 

208 

209 

210 
21 1 
212 

213 

214 

215 

216 

217 

218 

219 

220 

221 

nnc 

jC-.j~.uJ 

226 

'-\~7 
X- / 

223 

27.0 

231 


234 


235 

236 

237 

239 

240 

241 

245 

246 

247 
243 
249 


150 

251 


256 

260 


DIM X (301 , Y (30) , XG (360) ,YG(360> ,RDN(30,30> ,ZLEV(20> 
DIM R<6> , XNUM (12) , YNUM (12) ,NXD (30) ,NYD(30> 

DTRD-3 . 14159/180. 

DIM YM (70) ,Z (30,30) 

DAT A 0 , 30 , 60 , 90 , 1 20 , 1 50 ,180,210, 240 , 270 , 300 , 330 
FOR 1=1 TO 12 
NEXT I 


NLEV= 1 0 
N I NCR= . 05 


GOSUB 1000 
GOSUB 1600 

REM INITIALIZE F'LOTER 
0F’EN#9 , " . RS232" 

F'RINT#9; CHR£ (27) ; 132; ; 17: " : PR I NT #9 ; CHR-S ( 27 ) ; " 

PI X$=CONV-£ (PI X ) ; P2X$=C0NV-£ (F'2X ) : P1Y$=C0NV$ (PI Y) 
P2X=P1X+(X1-X0)*1010:P2Y=P1Y+ (Yl-YO) *1010 


N; 19; " ; “IN; " 
F'2Y^=C0NV-$ : (F'2Y ) 


I F X 0= X 1 THEN P 1 X =500 : P2 X = 1 0000 ; P 1 Y=600 : F'2 Y=7500 

F'l X-£=CONV£ (F'l X ) : F'2XT=C0NV* (F‘2X ) : P1Y4!=C0NV$ (P1Y) : F‘2Y-£=CQNV-£ (P2Y) 
F'FIELDT=P1X*+" , "+P1Y$+" , "+F'2X*+" , "+P2Y$ 

F'RINT#9; " IP" ; F'FIELDf ; " ; " 

GOSUB 220: RETURN 

REM SUBROUTINE TO SCALE AXISES 

REM X 0 , X 1 , Y 0 , Y 1 DEFINE STARTING AND ENDING POINTS IN INCHES 

REM XMI N , XMAX , YMIN , YMAX ARE THE MAX AND MIN VALUE ON AXISES 

REM SCALE X , SCALEY=SCALED DATA FOR TIC MARKS 

REM ND1VX ,NDIVY=NUMBER OF DIVISIONS ON AXISES 

DELTAX= (X 1— XO) /NDIVX 

DELTAY= ( Y 1-YO) /NDIVY 

SCALEX ( 1 ) =X0 : SCALEY ( 1 ) =Y0 

FOR ISUB=1 TO NDIVY: SCALEY ( ISUB+1 ) =SCALEY ( ISUB) +DELTAY: NEXT 
FOR I SUB= 1 TO NDIVX: SCALEX ( ISUB+1 )=SCALEX( ISUB) +DELTAX: NEXT 
SC4r=C0NV* (XO) +" , "+CONV.T ( X 1 ) +" , "+CONV* ( YO) +" , "+CONV* ( Y1 ) 

PR I NT#9 ; " SC " ; SC$ ; " ; " : RETURN 

REM SUBROUTINE TO SET SYMBOL MODE 

SYM 1 4:=C0NV^ ( W ) + " , " +CONV$ ( H ) 

F'RINT#9; “SI " ; SYM1*; "SM“ ; SYMf ; “ ; " 

RETURN 

REM THIS SUBROUTINE PLOTS SYMBOL SYM* AT LOCATION X,Y 

REM F'D=0 FOR PEN UP; PD=1 FOR PEN DOWN; LINET=LINE TYPE (0-4) 

REM L I NEL= LENGTH (PERCENT OF DIST. BETWEEN PI AND P2) 

IF PD=0 THEN PRINT #9; "PU; “ ; : ELSE PRINT#9; "PD; " ; 

IF PD=0 GOTO 239 
IF LINEL=0 THEN LINEL=4 
LT*=CONV* (LINED +" , “+CONV* (LINED 
F'R I NT #9 ; “LT“ ; LT* ; “ ; “ 

SYM2$=C0NV$ (X)+“ , “ +CONV.& ( Y ) :SYMB£=SYM$ 

PR I NT #9 ; “ PA “ ; SYM2* ; “ ; “ 

RETURN 


REM PLOT AXIS 


X 1 •■£=CONVJ : (XI) : Y1*=C0NV£ 
H PD;PA";X1$; " , Y0$; "; 

" ; "F‘D; PA" ; X0$; " , Yl$; 

" , " Y0£ ; " ; " 


( Y 1 ) 

" ; "PU; " 
" ; "PU; " 


F'RINT#9; "PU; " 

XO-£=CONV-$ (XO) : YO£=CONVT (YO) 

PRINT49; "PA" ; XO-S; " , " ; YO-S; " ; 

PRINT#9; "PA" ; X0$; " , "; Y0$; " ; 

F'R I NT #9 ; ” TL 1 . 0 , 0 . 0 ; PA " ; XO* ; 

FOR I SUB= 1 TO NDIVX+1 : SCALE$=CONVf (SCALEX ( ISUB) ) 

F‘R I NT #9 ; "F'A"; SCALED; " , " ; YO* ; "; XT; " : NEXT 
PRINT49; " TL 1 . 0,0.0; PA" ; XO*; " , "Y0$; " ; " 

FOR I SUB= 1 TO ND I VY+ 1 : SCALEr=CONV^( SCALEY ( ISUB) ) 
F'RINT#9; "F'A " ; XO-f ; SCALED; " ; YT ; " : NEXT 

RETURN 

REM SUBROUTINE TO PLOT NUMBER AXNUM AT LOCATION AX X , AYY 
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2o; , REM DIRUN, DIRISE: SET DIRECTION . 1,0 FOR HORIZONTAL 

263 REN SAW, SAH SET WIDTH AND HEIGHT IN CM. 

264 IF SAW=0 THEN SAW=0. 1 : SAH=0. 15 

266 AXX*=CONV* ( A X X ) : AYY*=CONV* ( AYY ) : AXNUM*=CONV* ( AXNUM) 

267 PRINTER; "SI" ; CONV* (SAW) ; " , " ; CONV* <SAH> ; " ; " 

268 PRINTER; "F'A" ; AX XT; " , " ; AYY* ; ” ; LB" ; AXNUM* ; CHR* (3) 

26R RETURN 

270 REM SUBROUTINE TO PLOT TEXT TXT* AT LOCATION TXX , TYY 

271 REM LOCATION IS THE LOWER LEFT HAND CORNER 

272 REM DIRUN, DIRISE; SET DIRECTION . 1,0 FOR HORIZONTAL 

273 REM SAW , SAH SET WIDTH AND HEIGHT IN CM. 

274 IF SAW=0 THEN SAW=0. 1 : SAH=0. 15 

275 PRINTER; "F'U; " 

276 TXX*=CONV* (TXX) : TYY*=CONV* (TYY) 

277 F'R I NT #9 ; "SI " ; CONV* ( SAW CONV* ( SAH ) ; " ; " 

273 PR I NT #9 ; "DI " ; CONV* (DI RUN) ; " , " ; CONV* (DIRISE) ; " ; " 

27R PRINTER; "PA" ; TXX*; " , " ; TYY*; " ; LB" ; TXT*; CHR* (3) 

280 RETURN 

281 REM SUBROUTINE TO CENTER DIGITS ON X AXIS 

282 NUMDIGIT=LEN (CONV* (AXNUM) ) 

283 AXX=AXX-(NUMDIGIT*SAW> /2/2. 54 

284 RETURN 

285 REM SUBROUTINE TO CENTER DIGITS ON Y AXIS 

286 AYY=AYY-SAH/2/2. 54; RETURN 

287 REM PLOT A CURVE OF NF' POINTS IN XG,YG 
283 XF-XG (0) : YF-YG (0) 

28R SYM2*=C0NV*(XP>+" , "+CONV*<YP> 

290 PR I NT #9 ; " F'U ; F'A " ; SYM2* ; " ; PD ; " 

2R1 FOR IF'L0T=1 TO NP 

292 XF'=XG ( I PLOT) 

293 YP=YG ( I PLOT) 

294 SYM2*=C0NV*(XF')+" , "+CONV*(YF') 

295 PRINT#9 ; "F'A" ; SYM2*; " ; F'D; " : NEXT I PLOT 

296 F'RINT#9; "F'U; " 

297 RETURN 

299 RETURN 

300 REM SUBROUTINE TO GET ERRORS ON HP 

301 F'RINT#9; CHR* (27) ; " . E; " 

302 , F'RINT#9; "OE; INPUT#9; ERROR*: PRINT ERROR* 

303 RETURN 

400 REM SUBROUTINE TO PUT AZIMUTH TIC MARKS 
410 TXT*="_" 

420 SAW=0. 1 : SAH=0. 1 
430 DIRUN=CSIN 

435 IF ABS (DIRUN) <0.01 THEN DIRUN=0. 

440 DIRISE=SINE 

445 IF ABS (DIRISE) TO. 01 THEN DIRISE=0 

450 TXX=XG(K) : TYY=YG <K ) 

460 GOSUB 270 
470 RETURN 

1000 REM INITIALIZE PLOTTER 

1010 P 1 X= 1 500 

1020 PI Y= 1250 

1030 X0=0 

1040 Y0=0 

1050 X 1=6 

1060 Y 1 =6 

1070 SAW=0. 1*2.54 

1080 SAH=0. 13*2. 54 

1090 CX=3. 0 

1100 CY=3. 0 

1110 XM I N=0 . O 

1120 XMAX=6.0 

1130 YMIN=0. 00 

1140 YMAX=6. 



,1 1*0 
1 1 70 
1 180 
1200 

1205 

1 206 
1210 
1220 
1 250 
1240 
1 250 
1255 
1260 
1270 
1280 
1 290 
1295 
1390 
1400 
1410 
1415 
1420 
1425 
1430 
1440 
1510 
1530 
1 600 
1610 
1620 
1 630 
1640 
1650 
1655 
1660 
1670 
1680 
1690 
1700 
1710 
1720 
1730 
1 800 
1810 
1820 
1830 
1840 
1850 
1360 
1870 
1880 
1910 
1920 
1925 
1930 
4900 
4910 
4920 

5000 

5001 

5002 

5003 

5004 

5005 
500 A 


Ul=i 

H=l 

GOSUB 225 A- 31 

GEM F'LDT BOUNDARY 

XNUM ( 1 ) =0: YNUM ( 1 ) =0 

DX=1 : DY= 1 : NP=5 

XG ( 1> =-2+CX : YG ( 1 ) =-2+CY 

XG <2) = <-2+CX : YG (2) =YG (1) 

X6 (3) =XG (2) : YG (3) =+2+CY 
XG (4) =XG < 1 ) : YG (4) =YG (3) 

X G < 5 ) = X G ( 1 ) : YG (5) =YG ( 1 ) 

PR I NT #9; "PU" 

FOR k= 1 TO NP 

XF’=XG <k) : YP=YG <K) 

SYM2$=C0NV$ (XP) +" , "+CONV* (YP) 

PRINT#9; "PA" ; SYM2$ ; ";PD;PA; " 

NEXT K 

SAW=. 25: SAH=. 25 
TXT ■£=" TAU = 60.0" 

TXX=2. 25: TYY=CY-2. 25 
D I RUN= 1 . 0 : D I R I SE=0 . 0 
GOSUB 270 
TXT*=" TAU = 0.0" 

TXX=2. 25: TYY=CY+2. 25 
GOSUB 270 

FOR P=1 TO 100: NEXT P 
RETURN 

REM INPUT INTERPOLATED BDN FIELD 

PRINT" INPUT FILE NAME OF INTERPOLATED DATA" 

INPUT D* 

0F'EN#2 , " / AREA 1 /CFOS . PLOT/ "+D$ 

1 = 1 

READ#2, 1 ; NX , NY , MIN , MAX 
PRINT NX , NY, MIN, MAX 
FOR J=1 TO NX 
FOR K=1 TO NY 
1 = 1 + 1 

READ#2 , I ; NXD ( J ) , NYD (K) ,RDN(J,K) 

PRINT J , K. , N X D ( J ) , NYD <K) , RDN < J , K ) 

NEXT K 
NEXT J 

REM SET UP CONTOUR ARRAYS 
ZMIN=MIN 
ZMAX=MAX 
FOR J=1 TO NX 
FOR F = 1 TO NY 
Z ( J , 1.) =RDN < J ,K) 

NEXT K 
NEXT J 

NXP=NX:NXP1=NX-1 
NYF-NY: NYP1=NY-1 
FOR N= 1 TO NLEV 

ZLEV ( N ) = I NT ( ZM I N* 1 00 ) / 1 00+N*N I NCR- 
NEXT N 
GOSUB 4900 
I DUB=0 

DELT AT=4/NYP 1 
DELTAR=4/NXP1 
FOR 1=1 TO NYP1 
T 0= ( 1-1 ) *DELT AT 
T 1 =T O+DELT AT 
FOR J=1 TO NXP 1 
R0= ( J— 1 ) *DELTAR 
R 1 =RO+DELT AR 



A- 32 


5003 1 

5009 

50 1 0 
50 1 1 

50 1 2 

5013 

5014 

5015 

50 1 6 

5017 

50 1 8 

50 1 9 

5020 

5021 

5022 

5023 

5024 

5025 

5026 

5027 
5023 

5029 

5030 

5031 

5032 

5033 

5034 

5035 


5036 


5037 
5030 
5082 
5033 
5035 
5086 
5090 
5 1 00 
5200 
5250 
5300 
5310 
5320 
5330 
5340 
5350 
5360 
5365 
5370 
5380 
5390 
5400 
5410 
5420 
5430 
5440 
5450 
5460 
5470 
54RO 


H3=Z < J + l ,1+1) 

H4 = Z ( J , 1 + 1 ) 

REM IF I = NYP 1 THEN H3=Z<J+1,1> 
REM IF I = NYP 1 THEN H4=Z(J,1) 
FDR h = 1 TO NLEV 
ZLEV=ZLEVO ) 

IF Hl-ZLEV GOTO 5023 
IF H2 = ZLEV GOTO 50 1 9 
IF H3.'=ZLEV GOTO 5020 
IF H4-ZLEV GOTO 5030 
GOTO 5100 

IF H3ZLEV GOTO 5021 -.ELSE GOTO 
IF H4'' ZLEV GOTO 5031 -.ELSE GOTO 
IF H4-.ZLEV GOTO 5033: ELSE GOTO 
IF H4ZLEV GOTO 5035:ELSE GOTO 
IF H2 = ZLEV GOTO 5026 
IF H3 : =ZLEV GOTO 5027 
IF H4< ZLEV GOTO 5036: ELSE GOTO 
IF H3< ZLEV GOTO 5028: ELSE GOTO 
IF H4YZLEV GOTO 5037: ELSE GOTO 
IF H4TZLEV GOTO 5032: ELSE GOTO 
IF H4 :• = ZLEV GOTO 5090 
TA=TO+ (Hl-ZLEV) / (H1-H4) *DELTAT 
) *DELTAR : GOTO 5080 


5022 

5032 

5034 

5036 


5035 

5029 

5033 

5031 

TB=T 1 : RA=RO: RB=RO+ < H4-ZLEV ) / (H4-H3 


TA=TO+ (H2-ZLEV) / (H2-H3) *DELTAT: TB=T1 : RA=R1 : RB=RO+ (H4-ZLEV) / (H4-H3 


) *DELTAR: GOTO 5080 

TA=TO+ ( H 1 -ZLEV) / < HI -H4 ) *DELTAT : TB=TO+ < H2-ZLEV ) / ( H2-H3 ) *DELTAT : RA= 


RO : RB=R 1 : GOTO 5080 

TA=T0:TB=T0+ (H2-ZLEV) / (H2-H3) *DELTAT : RA=RO+ < HI -ZLEV) / (H1-H2) *DELT 


AR : RB=R 1 : I DUB=0 : GOTO 5080 
I DUB= 1 : GOTO 5030 

TA=TO: TB=T 1 : RA=RO+ (Hl-ZLEV) / (H1-H2) *DELTAR: RB=RO+ (H4-ZLEV) / (H4-H3 
) *DELTAR: GOTO 5080 

TA=TO:TB=TO+ (Hl-ZLEV) / (H1-H4) *DELTAT: RA=RO+ (Hl-ZLEV) / (H1-H2) *DELT 
AR: RB=RO: IDUB=0: GOTO 5080 


I D LI E< = — 1 : GOTO 5031 

XA=RA+1 :‘XB=RB+1 : YA=TA+ 1 : YB=TB+1 


XA^=CONV^(XA) : XB$=CONV-$ ( XB) : YA$=CONV^ ( YA) : YBJ : =CONVf ( YB) 
PRINT#9; "PU" ; XA-T; " , "YA$; " ;PD";XB$; " ,"YB^; ":F'U; " 

IF IDUB=-1 THEN 5036 
IF IDUB= 1 THEN 5033 
NEXT K 
NEXT J 
NEXT I 
END 

FOR J=1 TO NXP 1 


FOR k=l TO NYP 1 

IF Z ( J , K ) = Z M I N THEN JMIN=J : KMIN=K 
IF Z ( J , K ) = Z M A X THEN JMAX=J : KMAX=K 
NEXT K 
NEXT J 

PRINT JMIN , KM I N , JMAX , KMAX 

D I RUN= 1 : D I R I SE=0 

RMIN= ( JMIN-1 ) *DELTAR 

RMAX= ( JMAX-1 ) *DELTAR 

TMIN= (KMIN-1 ) *DELTAT 

TMAX= (KMAX-1 ) *DELT AT 

XMIN=RMIN+1 

XMAX=RMAX+ 1 

YMIN=TMIN+1 

YMAX=TMAX+ 1 

TXX=XMIN: TYY=YMIN: TXTf="L" 

GOSUB 270 

TXX=XMAX : TYY=YMAX : TXT^="H" 
i;nci io n 7A 



5510 

5520 

5530 

5540 


TYY=~. 50 

TXTT="L = "+CGNV-S (ZMIN) +" H = " +C0NV* ( ZMAX ) + " INTERVAL = 
EV (2) -2LEV ( 1 ) ) 

GOSUB 270 ' 

END 


" +CONV£ ( ZL 
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) 


) 


) 

1 00 
1 10 
1 15 
120 
130 
200 
201 
202 

203 

204 

205 

206 
207 
203 
20 ? 
210 
21 1 

213 

214 

215 

216 

217 

218 

219 

220 

221 

One 
x- <_J 

226 

xL x. / 

228 

230 

231 



236 

237 

239 

240 

241 

245 

246 

247 

248 

249 

250 

" tcr . 1 
X_ v_l 1 

nc 

253 

254 

nee 
x~ U 

nc * 

wi£> 


DIM NA(30) , A2 <30) ,XG(360> ,YG(360) ,BDN<30,30) , ZLEV(20) 
DIM R ( 6 ) , X NUM (12), YNUM (12) 

DTRD=3. 14159/180. 

DIM AZ M < 30 ) , Z ( 30 , 30 ) 

DAT A 0 , 30 , 60 , 90 , 1 20 , 1 50 , 1 80 ,210, 240 , 270 , 300 , 330 
FOR 1=1 TO 12 
READ AZM(I) 

NEXT I 


NLEV=10 
N I NCR-. 50 


GOSUB 1000 
GOSUB 1600 

REM INITIALIZE F'LOTER 
OPEN# 1 ?, " . RS232 " 


F'RINT#9; CHR* (27) ; 132; ; 17: " : PRINT#9; CHR* (27) ; ".Ns 19; "; "IN; " 

P 1 X*=C0NV* ( P 1 X ) : P2 X* =CQNV* ( P2X > : P 1 Y*=C0NV* ( P 1 Y ) : P2Y*=C0NV* < F'2Y ) 


F'2X=P1X+ (X1-X0) >1010: P2Y=P1Y+ (Y1-Y0) *1010 


IF X0=X 1 THEN F’l X=500: F'2X = 10000: PI Y-600: F'2Y=7500 

P 1 X *=C0NV* ( P 1 X ) : F'2X*=C0NV* (F'2X ) : F'l Y* =C0NV* ( P 1 Y ) : P2Y*=C0NV* ( F'2Y ) 
F'FIELD*=P1X*+" , "+F'l Y*+" , "+P2X*+" , "+P2Y* 

PRINT#*?; "IP";F‘FIELD*; " ; » 

GOSUB 220: RETURN 


REM SUBROUTINE TO SCALE AXISES 

REM X 0 , X 1 , Y 0 , Y 1 DEFINE STARTING AND ENDING POINTS IN INCHES 

REM XMIN , XMAX , YMIN , YMAX ARE THE MAX AND MIN VALUE ON AXISES 

REM SCALEX , SCALEY=SCALED DATA FOR TIC MARKS 

REM NDIVX , NDI VY=NUMBER OF DIVISIONS ON AXISES 

DELTAX= ( X l-XO) /NDIVX 

DELTA Y= < Yl-YO) /NDI VY 

SCALEX ( 1 ) =X0: SCALEY ( 1 ) =Y0 

FOR I SUB= 1 TO NDI VY : SCALEY ( ISUB+1 ) =SCALEY ( ISUB) +DELTAY : NEXT 
FOR I SUB= 1 TO NDIVX: SCALEX ( ISUB+1 )=SCALEX< ISUB) +DELTAX: NEXT 
5C*=C0NV* (XO) +" , "+CONV* <X1 ) +" , "+CONV* (YO) + " , "+CQNV* <Y 1 ) 

PR I NT #9 ; " SC " ; SC* ; " ; " : RETURN 

REM SUBROUTINE TO SET SYMBOL MODE 

SYM1*=C0NV* (W) +" , " +CQNV* (H) 

PRINT#9; "SI"; SYM1*; "SM" ;SYM*; " ; " 

RETURN 

REM THIS SUBROUTINE PLOTS SYMBOL SYM* AT LOCATION X,Y 

REM F'D=0 FOR PEN UP; PD=1 FOR PEN DOWN; LINET=LINE TYPE (0-4) 

REM LINEL= LENGTH (PERCENT OF DIST. BETWEEN F'l AND F'2) 

IF F'D=0 THEN F'RINT#9; "F'U; " ; : ELSE PRINT#9; "PD; " ; 

IF PD=0 GOTO 239 
IF LINEL=0 THEN LINEL=4 
LT*=CONV* ( L I NET ) + " , " +CONV* ( L I NEL ) 

F'RINT#9; " LT " ; LT * ; " ; " 

SYM23 =CONV* < X ) + " , "+CONV* <Y> ; SYMB*=SYM* 

F'RINT#9; "PA" ; SYM2*; " ; " 

RETURN 

REM PLOT AXIS 


PRINT#9; "F'U; " 

XO*=CONV* (XO) : YO*=CONV* ( YO) : X 1 *=CONV* ( X 1 ) : Y1*=C0NV* (Y1 ) 


" ; "PD; PA" ; XI *; " , " ; YQ*; " ; " ; "F'U; " 

" ; "PD; F'A" ; XO*; " , ";Y1*; F'U; " 

" , " YO* ; " ; " 


F‘RINT#9; "PA" ; XO*; " , " ; YO*; " 

PR I NT #9 ; "PA" ; XO*; " , " ; YO*; " 

PR I NT #9 ; " TL 1 . 0 , 0 . 0 ; PA " ; XO* 

FOR I SUB= 1 TO NDIVX+1 : SCALE*=CONV* (SCALEX ( ISUB) ) 
PRINT#9; "PA" ; SCALE*; " , " ; YO*; " ; XT; " : NEXT 
PRINT49; "TL1 . 0 , 0. 0; PA" ; XO*; " , "YO*; " 

FOR I SUB= 1 TO ND I VY+1 : SCALE*=CONV* ( SCALEY ( I SUB ) ) 
PRINT#9; "F'A" ; XO* ; " , SCALE*; " ; YT; “ : NEXT 
RETURN 
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26 i • REM LOCATION IS THE LOWER LEFT HAND CORNER 

2e>2 REM DI RUN, DIRISE; SET DIRECTION . 1,0 FOR HORIZONTAL A-35 

263 REM SAW , SAH SET WIDTH AND HEIGHT IN CM. 

264 IF SAW=0 THEN SAW=0. 1 : SAH=0. 15 

266 AXX#=CONV# <AXX) ; AYY*=CONV# < AYY ) : AXNUM*=CONV# < AXNUM) 

267 PRINTER; "SI ” ; CONV# (SAW) ; ” , " ; CONV* (SAH) ; " ; " 

268 PRINT49; "PA" ; AXXT ; " , " ; AYY#; " ; LB" ; AXNUM*; CHR* (3) 

269 RETURN 

270 REM SUBROUTINE TO PLOT TEXT TXT* AT LOCATION TXX , TYY 

271 REM LOCATION IS THE LOWER LEFT HAND CORNER 

272 REM DIRUN, DIRISE; SET DIRECTION . 1,0 FOR HORIZONTAL 

273 REM SAW, SAH SET WIDTH AND HEIGHT IN CM. 

274 IF SAW=0 THEN SAW=0. 1 ; SAH=0. 15 

275 PR I NT 1*9; "PU; " 

276 TXX#=CONV# <TXX> : TYY*=CONV* <TYY) 

277 PR I NT#9 ; "SI " ; CONVT ( SAW CONV# < SAH ) ; " ; " 

278 • PRINT#9; "DI " ; CON Vi (DIRUN) ; " , " ; CONV* (DIRISE) ; "; " 

279 PRINT49; "PA" ; TXX*; " , ";TYY#; " ; LB" ; TXT#; CHR# (3) 

280 RETURN 

281 REM SUBROUTINE TO CENTER DIGITS ON X AXIS 

282 NUMDIGIT=LEN(CONV# (AXNUM) ) 

283 AXX=AXX- <NUMDIGIT*SAW> /2/2. 54 

284 RETURN 

285 REM SUBROUTINE TO CENTER DIGITS ON Y AXIS- 

286 AYY=AYY-SAH/2/2. 54: RETURN 

287 REM PLOT A CURVE OF NP POINTS IN XG , YG 

288 XP=XG<0> : YF'=YG(0) 

289 SYM2#=C0NV# (XP) +" , "+CONV# (YP) 

290 PR I NT #9 : "F’U; PA" ; SYM2#; " ;PD; " 

291 FOR IF'L0T=1 TO NP 

292 XP=XG ( I PLOT ) 

293 YF'=YG ( I PLOT ) 

294 SYM2#=C0NV# <XP) + " , "+CONV# (YP) 

295 F'RINT#9; "PA" ; SYM2#; " ; F‘D; " : NEXT IPLOT 

296 PRINT#9; "PU; " 

297 RETURN 

299 RETURN 

300 REM SUBROUTINE TO GET ERRORS ON HP 

301 F’R I NT #9 ; CHR# ( 27 ) ; ".E; " 

302 PR I NT#9 ; " OE ; " : I NPUT#9 ; ERROR# : PR I NT ERROR# 

303 RETURN 

400 REM SUBROUTINE TO PUT AZIMUTH TIC MARKS 

410 txt#="_" ORIGINAL PAGE BS 

420 SAW=0. 1 : SAH=0. 1 af POO^ QUAJJT02 

430 D I RUN=CS I N 

435 IF ABS (DIRUN) -0. 01 THEN D I RUN-0. 

440 DIRISE=SINE 

445 IF ABS < D I R I SE ) <0.01 THEN DIRISE=0 
450 T X X = X G ( K ) : T Y Y=YG ( K ) 

460 GOSUB 270 
470 RETURN 

1000 REM INITIALIZE PLOTTER 

1010 PI X= 1500 

1020 P 1 Y= 1 250 

1030 X0=0 

) 040 Y0=0 

1050 XI =6 

1060 Y 1 =6 

1070 SAW=0 . 1 *2. 54 

1080 SAH=0 . 13*2. 54 

1090 CX=3.0 

1100 CY=3.0 

1110 XMIN=0. 0 

1120 XMAX-6.0 


1 1 Ti 


VMTM-A 



A- 36 


1 1 3,0 
1 160 
1 170 
1130 
1200 

1 205 

1206 
1210 
1220 
1 230 
1 235 
1240 
1250 
1265 
1270 
1230 
1 290 
1 300 
1310 
1315 
1 330 
1340 
1350 
1360 
1390 
1400 
1410 
1415 
1420 
1425 
1430 
1440 
1450 
1455 
1460 
1470 
1480 
1490 
1500 
1510 
1520 
1530 
1 600 
1610 
1620 
1 630 
1640 
1650 
1655 
1 660 
1670 
1630 
1690 
1700 
1710 
1720 
1730 
1 300 
1810 
1820 
1830 
1840 
1 850 
1360 
1870 


GOSUb 2u'j 

W=1 

H= 1 

GOSUB 225 
REM PLOT BOUNDARY 
XNUM ( 1 ) -0: YNUM ( 1 ) =0 
DX= 1 : DY= 1 : NF’=360 
FOR 1=1 TO 6 
R ( I ) =1*0. 45 
FOR 1=0 TO 360 STEP 5 

CS I N=CQS ( 1 *DTRD ) : S I NE=S I N ( K*DTRD ) 

XG (1 )=R (I) *CSIN+CX 
YG (1 ) =R ( I ) *S I NE+CY 
SYM2-T=C0NV* <XP) +•• , "+CONV$(YP) 

XP=XG (k) : YP=YG(K) 

SYM2T=C0NVS ( XP ) +" , "+CONVS (YP) 

PR I NT #9 ; " PU ; PA " ; SYM2-S ; " ; PD ; " 

FOR P=1 TO 12 

IF K=AZM(P) THEN GOSUB 400 
NEXT P 
NEXT K 
NEXT I 
1=6 

FOR K=0 TO 330 STEP 30 
SAW=. 25: SAH=. 25 
TXTT= " "+CONV-J <K) 

DIRUN=COS(K*DTRD) 

IF AB3 (DIRUN) <0. 01 THEN DIRUN=0. 
DIRISE=SIN (k*DTRD) 

IF ABS (DIRISE) 70. 01 THEN DIRISE=0. 

TXX=R (6) *COS (1>DTRD) +CX 
TYY=R ( 6 ) *S I N ( K *DTRD ) +CY 
IF L.-:=90 THEN 1500 
IF K-=270 THEN 1500 
DIRISE=-DIRISE 
D I RUN=-D I RUN 

TXX=TXX+3*SAW*C0S <K*DTRD) 
TYY=TYY+3*5AH*SIN (K*DTRD) 

PRINT l:,TXX,TYY, DIRUN, DIRISE: GOSUB 270 
FOR P=1 TO 100: NEXT P 
NEXT F 
RETURN 

REM INPUT INTERPOLATED BDN FIELD 

PRINT" INPUT FILE NAME OF INTERPOLATED DATA" 

INPUT D$ 

0F’EN#2 , " / AREA1 / CFOS . PLOT / "+D£ 

1 = 1 

READ#2 , 1 ; NNAD , NAZM , FLUXUP , MIN , MAX 
PRINT NNAD , NAZM , FLUXUP , MIN , MAX 
FOR J=1 TO NNAD 
FOR K= 1 TO NAZM 
1 = 1 + 1 

READ#2 , I;NA(J) ,AZ <K) ,BDN(J,K) 

PRINT J , K , N A ( J ) , AZ <K> ,BDN(J,K) 

NEXT K 
NEXT J 

REM SET UP CONTOUR ARRAYS 
ZM I N=M I N 
ZMAX=MAX 
FOR J=1 TO NNAD 
FOR 1=1 TO NAZM 
Z (J ,K) =BDN ( J ,K) 

NEXT K 
NEXT J 

NRAD=NMAD : NRM 1 =NRAD— 1 



t 91 (') 
1920 
1925 
1 920 
4900 
4910 
4920 

5000 

5001 

5002 

5003 

5004 

5005 

5006 

5007 

5008 

5009 

50 1 0 
501 1 

50 1 2 

5013 

5014 

50 1 5 

50 1 6 

5017 

5018 

5019 

5020 

5021 

5022 

5023 

5024 

5025 

5026 

5027 

5028 

5029 

5030 

503 1 

5032 

5033 


5034 

5035 


5036 


5037 

5080 


5082 

5083 

5085 

5086 
5090 
5 1 00 
5200 
5250 
5300 
53 1 0 
5320 
5330 


FOR N=1 ID NLEV 

ZLEV (N) =TNT (ZMIN*10) / 1 O+N+N I NCR 
NEXT N 
GOSUB 4900 
I DUB=0 

DELTAT=3. 14159/ (NTH 1 ) *2 
DELT AR=R ( 6 ) / ( NRT1 1 ) 

FOR 1=1 TO NTM1 
T0= (I-l) *DELT AT 
T 1 =T O+DELT AT 
FOR J = 1 TO NRM1 
R0= (J-l ) *DELT AR 
R 1 =RO+DELT AR 


H1 = Z (J , 1 ) 

H2=Z (J+l , I) 

H3=Z (J+l , 1+1 ) 

H4=Z (J , 1+1 ) 

IF I=NTM1 THEN H3=Z(J+1,1) 
IF I=NTM1 THEN H4=Z(J,1> 

FOR F=1 TO NLEV 
ZLEV=ZLEV (K) 

IF Hl;=ZLEV GOTO 5023 
IF H2.=ZLEV GOTO 5019 
IF H3' =ZLEV GOTO 5020 
IF H4ZZLEV GOTO 5030 
GOTO 5100 

IF H3-. ZLEV GOTO 5021: ELSE 
IF H4-' ZLEV GOTO 5031: ELSE 
IF H4 =-■ ZLEV GOTO 5033: ELSE 
IF H4 ZLEV GOTO 5035: ELSE 
IF H2 - =ZLEV GOTO 5026 
IF H3-=ZLEV GOTO 5027 
IF H4-- ZLEV GOTO 5036: ELSE 
IF H3'‘ ZLEV GOTO 5028: ELSE 


GOTO 5022 
GOTO 5032 
GOTO 5034 
GOTO 5036 


GOTO 5035 
GOTO 5029 


IF H4 ^ ZLEV GOTO 5037: ELSE GOTO 5033 
IF H4' y ZLEV GOTO 5032: ELSE GOTO 5031 
IF H4 * = ZLEV GOTO 5090 

TA=T0+ (Hl-ZLEV) / ( H 1 -H4 ) *DELT AT : TB=T 1 : R A=R0 : RB=R0+ (H4-ZLEV) / (H4-H3 
) *DELT AR : GOTO 5080 

TA=T0+ (H2-ZLEV) / < H2-H3 ) *DELT AT : TB=T 1 : RA=R 1 : RB=R0+ (H4-ZLEV) / (H4-H3 
) +DELTAR: GOTO 5080 


TA=TO+ (Hl-ZLEV) / ( H 1 -H4 ) *DELT AT : TB=T 0+ (H2-ZLEV) / (H2-H3) *DELTAT: RA= 
RO : RB=R 1 : GOTO 5080 


TA=TO: TB=TO+ (H2-ZLEV) / (H2-H3) *DELTAT : RA=RO+ (Hl-ZLEV) / (H1-H2) *DELT 
AR: RB=R1 : IDUB=0: GOTO 5080 
IDUB=1 : GOTO 5030 


TA=T0: TB=T 1 : RA=RO+ (Hl-ZLEV) / (H1-H2) *DELTAR: RB=RO+ (H4-ZLEV) / (H4-H3 
) *DELTAR: GOTO 5080 

T A=T 0 : TB=T 0+ (Hl-ZLEV) / (H1-H4) *DELTAT : RA=R0+ (Hl-ZLEV) / (H1-H2) *DELT 
AR: RB=RO: IDUB = 0:G0T0 5080 


IDUB=— 1 : GOTO 5031 

X A=RA+COS (TA) +CX : XB=RB*COS (TB) +CX : YA=RA*SIN <TA) +CY: YB=RB*SIN (TB) + 


CY 

XAT=CONVT (XA) 


XBT=C0NVT(XB> : YAT=CONVf ( YA) : YB$=CONV$ ( YB) 


PR I NT#9 ; "PU" ; X AT ; " , " YA?; 1 
IF IBLIB=-1 THEN 5036 
IF I DUB= 1 THEN 5033 
NEXT ) 

NEXT J 
NEXT I 
END 

FOR J=1 TO NNAD 
FOR K— 1 TO NAZN 

IF Z(J,K)=ZMIN THEN JMIN=J; 
IF Z(J,F1=ZMAX THEN JMAX=J: 


PD" ; XBT ; " , "YB?; ";PU; 


FNIN=K 

KNAX=F 
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5350 NEAT J 

5360 ' PRINT JMI N , KM I N , JMAX , !■ MAX 
5365 DIRUN=1 : DIRISE=0 
5370 RMIN= (JMIN-1 ) *DELTAR 
5380 RMAX= (JMAX- I ) *DELTAR 
5390 TMIN= O' MIN-1 ) *DELTAT 
5400 TMAX= <K MAX-1 ) *DELTAT 
5410 XMIN=RNIN*C0S<TMIN)+CX 
5420 XMAX=RMAX*COS (TMAX) +CX 
5430 YMIN=RMIN*SIN (TMIN) +CY 
5440 YMAX=RMAX*SIN (TMAX) +CY 
5450 TXX = XMIN: TYY=YMIN: TXTJ : ="L" 

5460 GOSUB 270 

5470 TXX = XMAX : TYY=YMAX : TXT$= a "H" 

5480 GOSUB 270 

5490 END 

5500 TXX=0.30 

5510 TY Y=— . 50 

5520 TXT$="L = "+CONV-P (ZMIN) + " H = "+CONV$ < ZMAX ) +" INTERVAL 
EV <2) -ZLEV ( 1 ) ) 

5530 GOSUB 270 

5540 END 


"+CONV* (ZL 



A- 39 


10 REM PROGRAM CFOS 

20 PRINT" CFOS BI-DIRECTIONAL REFLECTANCE INTERP. PROGRAM" 

45 DIM BDR (13,15) , BDN (30 , 30) , NAN (30) , AZ IN (30) 

50 DIM NAD ( 13, 15) , AZI ( 13, 15) 

60 DATA 2. 5, 7. 5, 12.5, 17.5,22.5,27.5,32.5,37.5,42.5,47.5,52.5,57.5 

70 DATA 15,45,75,105,135,165,195,225,255,285,315,345 

75 PI=3. 14159; P0VR2=PI/2 

100 FOR J=1 TO 30 

110 FOR )=1 TO 30 

120 BDN ( J , K ) =0 . 

150 NEXT K 

160 NEXT J 

165 PRINT "ENTER THE NUMBER OF NAD INTERPOLTION POINTS" 

167 INPUT NNAD 

170 FOR J=1 TO NNAD 

190 NAN (J) =PI /3/NNAD* (J-l ) 

200 NEXT J 

205 DTRD=3. 14159/180. 

210 DEF FN ACOS ( X ) =ATN (SQR (1-X*X)/X) 

215 DEF FN ACOS ( X ) =ATN <SQR ( 1-X*X > XX ) 

216 PRINT "ENTER THE NUMBER OF AZIMUTHAL INTERPOLATION POINTS" 

217 INPUT NAZM 

220 FOR J=1 TO NAZM 

230 AZ IN ( J ) =2*F' I / NAZM* ( J - 1 ) 

240 NEXT J 

250 PRINT" INPUT DATA FILE NAME " 

260 INPUT D-f 

270 0F'EN#3 AS INPUT , "/AREA 1 /CFOS. BDR/ "+D$ 

280 READ#3, 1 ; NS , FLUXUP , M IN , MAX 

290 1=1 

300 FOR J=1 TO NS 

310 FOR 1=1 TO 15 

320 1=1+1 

330 READ#3 , I ; J , h , NAD ( J , K) ,AZI (J,K) ,BDR(J,K) 

340 PRINT J , K , NAD ( J , K ) , AZ I < J , K ) , BDR < J , K ) 

350 NEXT K 

360 NEXT J 

370 PRINT" UPWARD FLUX DENSITY EQUALS ", FLUXUP 

380 PRINT" MIN AND MAX VALUES OF BDR FUNCTION ", MIN, MAX 

400 GOSUB 1000 

500 REM SUBROUTINE TO CALCULATE WEIGHT 
510 DELTA=PI/6. 

520 IF THETA y DELTA GOTO 570 

530 WGHT = 1 . 0E+06 

540 IF THETA- 1 0. 001 THEN RETURN 

550 WGHT=DELTA/ THETA 

560 RETURN 

570 WGHT=0. 

580 IF (THETA-DELTA) ;=DELTA THEN RETURN 
590 WGHT=2. O-THETA/DELTA 

600 RETURN 

700 REM SUBROUTINE TO INTERPOLATE IN BDR ARRAY 
710 SUMF=0. 

720 SUMW=0. 

730 FOR M=JMIN TO JMAX 
740 FOR N=KMIN TO KM AX 

750 Z I =NAD (M , N) 

760 A I =AZ I (M , N) 

770 THET A=COS ( Z I ) *COS ( Z PRM ) +SIN ( Z I ) *SIN ( ZPRM ) *COS ( A I -AF'RM ) 

780 IF THETAKO. GOTO 840 

790 THETA= FN ACOS (THETA) 
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8 U U 

’si 6 ’ 
820 
830 
840 
850 
860 
870 
880 
900 
1 000 
1010 
1020 
1030 
1040 
1050 
1060 
1070 
1075 
1080 
1085 
1090 
1 1 00 
1 1 20 
1130 
1 1 50 
1160 
1 170 
1 1 SO 
1 190 
1200 
1210 
1230 
1240 
1250 
1260 
1270 
1280 
1290 
1 300 
1310 
1320 
1330 
1340 
1 350 
1 360 
1400 
1410 
1420 
1430 
1440 
1450 
1460 
1470 
1480 
1485 
1490 
1495 
1500 
1510 
1520 
1530 
1 535 
1540 
1550 


GuSUB 500 

W=WGHT 
SUMW=SUMW+W 
SUMF=SUMF+W*BDR(M,N) 

NEXT N 
NEXT M 
INTERF-O. 

IF SUMW' =0. THEN RETURN 
I NTERF‘=SUMF / SUMW 
RETURN 

REM CALCULATE BDR VALUES AT REGULAR ARRAY POINTS 
FOR J=1 TO NNAD 
FOR K= 1 TO NAZM 
ZPRM=NAN < J ) 

APRM=AZ IN (K) 

REM FIND JMIN , JMAX , KM IN , KMAX TO LIMIT REGION OF INTERPOLATION 
S I ND=S I N < ZPRM ) *8 1 N ( P I -AF’RM ) 

COSD=SDR ( 1-SIND*5IND) *SGN (COS ( APRM) ) 

COSA=COS ( ZPRM ) /COSD 
IF COS A. =1.0 THEN C0SA=1.0 
IF COSAX =-1.0 THEN C0SA=-1 . 0 
A= FN ACOS (COS A) 

D=ATN(SIND/COSD*SGN (COS (APRM) ) ) +P0VR2 

AINC=INT(A/PI*lB>+7 

DINC=INT (D/F’I*18)-1 

JMIN=AINC— 1 

JMAX=A I NC+ 1 

1 MIN=DINC— 1 

KMAX=DINC+1 

IF JMIN;=1 THEN 1230 

JMTN=JMIN+1 

GOTO 1190 

IF JMAX\=13 THEN 1260 
JMAX = JMAX- 1 
GOTO 1230 

IF KM IN ' = 1 THEN 1290 
KM IN=KM IN+ 1 
GOTO 1260 

IF KMAX-: = 15 THEN 1320 

KMAX=KMAX-1 

GOTO 1290 

GOSUB 700 

BDN ( J , K ) = I NTERP 

PRINT BDN < J ,K) , JMIN , JMAX , KMIN , KMAX 
NEXT K 
NEXT J 

0F'EN#4 AS OUTPUT, " /AREA 1 /CFOS. PLOT/ "+D# 

MIN=99999 

MAX=0 

FOR J = t TO NNAD 
FOR K=1 TO NAZM 

IF BDN < J ,K) < MIN THEN MIN=BDN(J,K) 

IF BDN < J ,K) >MAX THEN MAX=BDN < J ,K) 

NEXT K 
NEXT J 

PRINT "MIN AND MAX OF INTERPOLATED FIELD ", MIN, MAX 
WR I TE44 , 1 ; NNAD , NAZM , FLUXUP , M IN , MAX 
1 = 1 

FOR J=1 TO NNAD 
FOR 1=1 TO NAZM 
1 = 1 + 1 

WRITE44, I ; NAN ( J ) ,AZIN(K) ,BDN(J ,K> 

PRINT NAN (J ) , AZIN <K) , BDN (J ,K ) 

NEXT h 
NEXT J 



END 
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10 REM CFOS BI-DIRECTIONAL REFLECTANCE STANDARD ERROR PROGRAM 
20 DIM N (30, 12, 12) ,NAN( 12) ,A2IM( 12) 

30 DIM NBAR ( 12, 12) , NSTD ( 12, 12) , SSN ( 12, 12) , SNS ( 12, 12) 

40 DIM DT (301 

100 PRINT" INPUT THE NUMBER OF INTERPOLATED DATA FILES TO USE" 

110 INPUT NUM 

130 FOR 1=1 TO NUM 

140 PRINT "INPUT THE DATASET NAME FOR FILE NUMBER", I 

150 INPUT • D£ ( I ) 

155 D$ ( I ) ="/AREAl /CFOS. PLOT/ "+D$ ( I ) 

160 NEXT I 

200 FOR 1=1 TO NUM 

210 0PEN#2 AS INPUT , D$ ( I ) 

220 10UNT=1 

230 FOR 3=1 TO 12 

240 FOR H =1 TO 12 

250 FOUNT =KOUNT + 1 

255 READ42, 1 ; NNAD , NAZM , FLUXUP , MIN , MAX 

260 RE AD#2 , FOUNT ; NAN ( J ) ,AZIM(K) ,N(I ,J,K> 

270 NEXT k 

2S0 NEXT J 

290 CL0SES2 

300 NEXT I 

310 FOR 3=1 TO 12 

320 FOR k=l TO 12 

330 SSN ( 3 , k ) =0 

340 SNS (3 , K) =0 

350 NBAR (3 , K) =0 

360 NSTD (3 , k ) =0 

370 NEXT K 

3S0 NEXT 3 

400 FOR 3=1 TO 12 

410 FOR 1=1 TO 12 

420 FOR 1=1 TO NUM 

430 SSN <3,k) =SSN<J,K)+N< I ,J,K>*N (I ,3,K) 

440 SNS (3 ,K) =SNS ( J , K ) +N ( I , J , K ) 

450 NBAR ( 3 , Y ) =NBAR < J , K ) +N < I , J , K ) /NUM 

460 NEXT I 

470 SNS ( 3 , K ) =SNS < 3 , K ) *SNS ( 3 , K ) 

480 NSTD (3 , K ) =SQR ( (SSN ( J , K) -SNS (3 , K) /NUM) / (NUM-l ) ) 

490 NSTD <3 ,K) =NSTD ( J ,K) /NBAR (3 , K) 

500 PRINT J , K , NSTD < J , K ) * 1 00 

510 NEXT K 

520 NEXT 3 

530 AVGERR=0 

540 FOR 3=1 TO 12 

550 FOR 1=1 TO 12 

560 AVGERR=AVGERR+NSTD < J , K) *100/144 

570 NEXT V 

580 NEXT 3 

590 PR I NT "OVERALL STANDARD ERROR FOR BDR" , AVGERR 
600 END 



